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ABSTRACT 


SDI  tachnologlas  may  affect  the  concepts  of  weapons 
systems  and  reshape  the  battlefields  of  the  future,  in  this 
research  a  few  of  the  SDI  systems  were  analyzed  and  modeled 
mathematically.  The  different  models  were  gathered  in  a 
software  package  that  may  be  run  on  a  personal  computer. 
The  intention  was  to  produce  a  handy  tool  for  preliminary 
studies  and  designs. 

The  following  topics  are  covered:  aerodynamic  design  of 
hyperveiocity  projectiles,  flight  simulation  of 
hypervelocity  projectiles,  railgun  simulation,  ground  track 
of  orbits,  guidance  and  homing. 

A  few  case  studies  were  analyzed  to  suggest  and 
demonstrate  potential  use  of  these  models. 
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I.  I£ 


In  March  23,  1983  tha  Praaidant  of  tha  Uni  tad  Stataa 
announcad  ha  waa: 

diracting  a  conprahanaiva  and  intanaiva  affort  to  dafina 
a  long-tam  raaaarch  and  davalopaant  prograa  to  achiava 
our  ultimate  goal  of  eliminating  tha  threat  poaad  by 
strategic  nuclear  miaailaa. 

In  1984  tha  Oapartmant  of  Dafanaa  aatabliahad  an 
organization  to  expand  and  accalarata  raaaarch  in  this  area. 
The  research  program  was  called  the  "Strategic  Defense 
Initiative,"  better  known  as  the  SDI. 

Since  that  day,  research  and  studies  on  tha  subject 
continue  at  an  accelerated  rata.  There  Is  also  a  growing 
controversy  over  the  potential  of  various  SDI  technologies 
and  tlie  possibilities  for  applying  them  in  affordable  weapon 
systems  that  would  be  effected  when  needed. 

Despite  the  controversy  over  the  SDI  system's  ability  to 
achieve  Its  strategic  role,  there  Is  a  wide  acceptance  that 
SDI  technologies  may  also  affect  the  concepts  of  weapons 
systems  and  reshape  tha  battlefields  of  the  next  century. 
It  is  therefore  Important  to  study  the  additional  non-SDI 
contributions  that  SDI  technologies  could  make  on  the 
future . 

In  this  research  a  few  of  the  SDI  weapon  systems  were 
analyzed  and  modeled  mathematically.  The  various  models 
were  gathered  In  a  software  package  that  may  be  run  on  an 
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in«xp«nalv«  dask-»ountad  parsonal  cof^putar.  Thia  anal/tioal 
tool  was  built  to  assura  an  aasy  aooaas  systaa.  It  will 
allow  tha  usair  to  laam»  to  study  and  avan  to  parfors  a 
praliainary  dasign  of  futura  %faapon  systaaitt  basad  on  thasa 
taohnologias.  Tha  aaphasis  was  to  provida  a  systaa  that 
will  allow  a  vary  quick,  intaraotiva  and  aasy-to-oparata 
working  tool. 

Most  nodals  wara  mrittan  in  ST-Basic.  Thay  maka  usa  of 
tha  rich  graphics  capabilities  of  tha  ATARI  ST-1040 
conputar.  Tha  nachina-languaga  version  was  conpilad  for 
■oat  nodals.  A  Fortran  version  of  tha  orbit  nodal  was 
produced  and  run  on  tha  IBM  aainfrana.  This  version  is  used 
whenever  very  long  cases  are  run,  and  tha  graphic 
presentation  is  not  needed. 

The  systea  was  built  in  nodules.  By  using  this 
approach,  it  becane  easier  to  add  nore  modules  or  to  improve 
the  existing  ones,  easily.  In  the  process  of  this  research, 
models  were  designed  on  the  following  topics: 

a.  Aerodynamic  design  of  hyper-velocity  projectiles. 

b.  Flight  simulation  of  hyper-velocity  projectiles. 

c.  Evaluation  of  the  railgun  performance. 

d.  Ground  track  of  orbits. 

e.  Guidance  and  homing. 

Other  models  for  study  of  these  topics  do  exist. 
However,  most  of  them  are  very  sophisticated  and  complex. 
The  majority  of  them  were  written  to  be  operated  on  large 
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MinfraMft.  Soi«  of  thm  oon  bo  uood  only  by  oxports  in 
thoso  spool fio  fioldo. 

Tho  opp;?oooh  of  this  work  wos  to  orooto  o  vory  sinplo 
and  handy  tool  to  uso.  Xt  was  proparod  for  uso  by  systoa 
onginoora  and  woapona  dosignors  that  nood  to  know  tho  "big 
pioturo"  without  going  dooply  into  tho  fino  dotails  of  oach 
cosponont.  This  approach  is  vory  isportant  in  tho 
proliainary  dosign  phaso,  in  which  aany  itoratlons  aro  dona 
until  a  woll-balancod  proliainary  dosign  is  aohiovod. 

Whonovor  a  aoro  accurato  and  rofinod  study  is  noodod,  a 
aoro  sophist icatod  aodol  has  to  bo  usod.  Tho  uso  of  this 
siaplo  aodol  will  assuro  tho  officiont  uso  of  tho 
sophisticatod  aodols  bocauso  tho  tasks  will  bo  bottor 
dofinod  and  tho  usor  will  havo  provious  oxporioneo  with  tho 
topics. 

In  tho  following  chaptors  oach  aodulo  of  the 
mathoaatical  aodols  will  bo  oxplainod.  A  fow  case  studies 
will  bo  given  in  order  to  suggest  and  doaonstrato  potential 
uso  of  thoso  aodols.  Conclusions  and  roconaondatlona  will 
close  this  research  report. 
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II.  THE  DRAG  MODEL  rAgRQDYNAMIC  COirEICI 


gaXIMMIQW) 


X.  GENERAL 

Th«  purpose  of  th«  Aorodynamic  Coaffieiants  Eatiaation 
aodal  is  to  provlda  initial  valuaa  for  tha  praliainary 
daaign.  Tha  dasignar  givas  as  an  input  tha  valua  of  tha 
main  gaoaatrical  paraaatars  that  dataraina  tha  shapa  of  tha 
body*  Tha  aodal,  basad  on  doubla  intarpolation  of 
aarodynaaic  graphs  produce  as  an  output  tha  valuas  ralavant 
to  tha  chosan  shapa.  Tha  prograa  also  providas  a  soalad 
schaaa  of  tha  body  and  tha  liaitations  of  tha  cantar  of 
gravity  location  to  aasura  stability. 

Aftar  studying  tha  output  tha  dasignar  aay  aithar  daeida 
to  changa  tha  input  paraaatars  and  to  start  a  nav  design 
iteration,  or  if  ha  is  satisfied  with  tha  results,  ha  aay 
proceed  to  the  next  phase  of  tha  design. 

B.  THE  MATHEMATICAL  MODEL 

1.  List  of  constants 

The  list  of  constants  used  in  this  aodal  is  given  in 
Table  2.1.  The  following  description  is  given  for  each 
constant: 

a)  Tha  text  symbol 

b)  The  definition 

c)  The  value 

d)  The  units. 
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TABLE  2.1 


TABLE  2.2 
LIST  OF  VARIABLES 

SYMBOLS  MEANINGS  UNITS 


Computer 

Code 

Text 

XI 

Xl 

X  coordinate  of  point 

1 

Non 

Dimen. 

X2 

X2 

X  coordinate  of  point 

2 

Non 

Dimen . 

X3 

X3 

X  coordinate  of  point 

3 

Non 

Dimen . 

YAl 

^A1 

Y  coordinate  of  point 
Curve  A 

1, 

Non 

Dimen . 

YA2 

Ya2 

Y  coordinate  of  point 
Cui-ve  A 

2, 

Non 

Dimen . 

12 
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»  ■  ■  j 

T^LE  2.2  (COMTINUED) 

SYMBOLS 

MEANINOS 

UNITS 

<• 

Computer 

Code 

Text 

j 

YA3 

^A3 

Y  coordinate  of  point  3, 
Curve  A 

1 

Non  Dimen.  ! 

j 

i 

YBl 

’^Bl 

Y  coordinate  of  point  i. 
Curve  B 

1 

Non  Dimen. 

i 

YB2 

Yb2 

Y  coordinate  of  point  2, 
Curve  B 

Non  Dimen. 

YB3 

Yb3 

Y  coordinate  of  point  3, 
Curve  B 

Non  Oimen. 

1 

YCl 

^Cl 

Y  coordinate  of  point 
Curve  C 

Non  Dimen. 

YC2 

Yc2 

Y  coordinate  of  point  2, 
Curve  c 

Non  Dimen. 

• 

YC3 

Yc3 

Y  coordinate  of  point  3, 
Curve  c 

Non  Dimen. 

YDl 

^dl 

Y  coordinate  of  point  1, 
Curve  l> 

Non  Dimen. 

YD2 

Yd2 

Y  coordinate  of  point  2, 
Curve  D 

Non  Dimen. 

YD3 

^d3 

Y  coordinate  of  point  3, 
Curve  D 

Non  Dimen. 

CD 

Cd 

Coefficient  of  drag 

Non 

AD 

d 

Cone  half  angle 

Deg 

AR 

e 

Cone  half  angle 

Rad 

RNRB 

Rr 

Ratio  of  radii 

Non  Dimen. 

L 

L 

The  cone  length 

m 

-  ■ 

RP 

R 

The  radius  of  cone  base 

m 

XCP 

Xcp 

Location  of  c.g. 

(measured  from  nose  tip) 

m 
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TABLE  2.2  (CONTINUED) 


SYMBOLS  MEANINGS  UNITS 


Computer 

Code 

Text 

VP 

Vp 

The  coii«<  volume 

in3 

MP 

Mp 

The  cone  mass 

)cg 

XS 

Xs 

X  coordinate  selected 

Non 

■  1  ;■ 

Dimen. 

YS 

Xa 

Y  coordinate  selected 

Non 

b'imen. 

ZS 

Zs 

Z  value  at  Xg,  Yg 

Ncn 

Dimen. 

X 

X  coordinate 

Non 

Dimen. 

Y 

Y  coordinate 

Non 

bimen. 

Z 

The  function  value 

Non 

Dlmen. 

3.  eagjg,  AggvHaptions  and  aelationg 

The  construction  of  this  model  is  based  upon  the 
following  basic  assumptions  and  relations: 

a.  For  this  preliminary  design,  only  the  particle 
dynamics  of  the  projectile  are  relevant.  In  reality 
vehicular  angular  motions  exist,  and  alter  the  new 
flow  field,  and  hence  also  affect  the  heat  transfer 
mechanism.  In  this  model  it  was  assumed  that  a  single 
average  value  of  the  coefficient  of  drag  may  properly 
represent  the  aerodynamic  characteristics  of  a  ncn- 
maneuvering  projectile  during  launch  and  reentry,  and 
the  earth's  atmosphere. 

b.  Since  most  of  the  cases  of  interest  are  related  to 
hypersonic  flight,  the  value  of  the  drag  coefficient 
is  Independent  of  the  Mach  number. 

c.  It  is  also  assumed  that  most  relevant  reentry  bodies 
may  approximate  a  blunt  cone.  (Atmospheric  flights  of 
hypersonic  speeds  demand  that  frictional  drag  be 
smaller  relative  to  the  pressure  drag  to  prevent 
overheating.)  The  engineering  implications  of  this 
demand  are  that  reentry  and  high  hypersonic  launched 
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projectiles  must  have  large  pressure  drag,  i.a.,  blunt 
shape . 

d.  The  stiape  of  blunt-cone  (see  Figure  l.l)  may  be 
satii^factorlly  defined  by  the  following  two 

'  parameters : 

a.  The  cone  half  angle  (6) 

b.  The  ratio  of  cone  nose  radius  to  base  radius 
(Rr)  • 

e.  The  function  of  the  drag  coefficient  vs.  cone 
bluntness  (for  a  constant  value  of  cone  half  angle) 
may  be  approximated  by  polynomials  ( fourth  degree 
polynomials  were  used  in  this  version) . 

f.  The  graph  of  cone  drag  vs.  cone  bluntness  of  [Ref.  1] 
was  chosen  to  represent  the  data  base  for  this  model. 

g.  The  function  of  the  drag  coefficient  (C^)  vs.  the  cone 
half  angle  (for  any  given  value  of  cone  nose  radius  to 
base  radius  ratio)  may  be  approximated  by  a 
polynomial. 

h.  The  calculation  of  the  location  of  the  center  of 
pressure  is  based  on  [Ref.  2]  where  an  analysis  was 
described  for  conical  noses.  Based  on  this  analysis 
it  is  assumed  that  the  location  of  the  center  of 
pressure  is: 


Xcp  -  (L) (2/3)sec2e 


where: 

Xqp  is  the  location  of  the  center  of  pressure  measured 
from  the  nose  tip 

L  is  overall  length  of  the  nose 

0  is  the  cone,  half  angle. 

This  is  considered  a  conservative  approach  that 
includes  a  safety  factor,  whenever  the  projectile 
includes,  besides  the  nose,  also  tall  and/or  fins. 

i.  An  estimate  of  the  projectile  volume  and  mass  is 
given.  The  estimate  is  based  on  the  volume  of  a  blunt 
cone.  The  mass  estimate  is  based  on  a  typical  value 
fcr  a  specific  weight  for  guided  warheads  [Ref.  3]. 
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Th*  data  base  for  this  program  was  provided  as  a  4 
curve-graph  (see  Figure  2.1).  Each  curve  represents  the 
function  of  the  drag  coefficient  (Cjj)  vs.  the  ratio  of  the 
cone  nose  radius  to  the  base  radius  for  a  Known  value  of 
cone  half  angle  (6).  These  curves  are  typical  of  other 
aerodynamic  carves  that  represent  the  value  of  an 
aerodynamic  variable  as  a  function  of  one  geometric  vaT  j.<able 
while  other  geometric  x'alues  are  held  constant.  This  is  the 
common  method  of  presenting  the  empirical  results  of 
aerodynamic  tests  that  normally  are  done  in  a  similar  method 
(rv\nning  over  values  of  one  parameter  while  holding  the 
value  of  other  parameters  constant) .  In  more  general  terms, 
it  is  a  case  of  function  f  that  relates  the  value  of  a 
variable  a  to  the  values  of  two  others,  x,y  (see  Figures  2.2 
and  2.3).  The  values  of  z  =  z(x)  are  given  for  a  few  values 
of  constant  value  of  y  (y  serving  as  a  parameter) . 

The  aerodynamic  data  in  this  specific  case  consists 
of  3  curves:  Curves  A,  B,  and  C.  Every  curve  is 
j^epresented  by  a  set  of  values  of  the  drag  coefficient  (z) 
versus  the  cone  angle  (y)  and  the  radius  ratio  (x)  .  For 
each  point; 

^Mn  **  f(%n#^Mn) 
where  the  notations  are: 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


RE-ENTRY  BODY  PARTICLE  MECHANICS 


.  I 


Figure  2.1  vs  Cone  Bluntness 
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N  •  th«  curv«  (A,  B,  or  C) 
n  »  Sorial  nuabor  of  tho  point. 

Th*  first  noduls  of  this  nodal  is  rssponsibls  to 
producs  for  sach  function  (raprassntsd  by  values  of  z  vs. 
values  of  x)  an  interpolating  polynonlal.  The  interpolating 
polynomials  were  given  in  the  Lagranglan  form.  (The 
Lagrangian  form  was  selected  because  it  is  suitable  for  the 
general  nonuni formly-spaced  x  values  case) . 

When  a  specific  value  of  Zg  ■  for  a  given 
value  of  X  >  Xg  and  y  «  is  needed,  the  next  module 
(first-interpd)  produces  first  the  functions  of  z  f(y) 
(see  Figure  2.3)  tor  the  selected  value  of  x  «  Xg.  This  is 
done  by  interpolation  of  all  the  Z  «  Z(x)  curves  at  the 
value  of  X  »  Xg. 

The  next  module  (second-interpd)  takes  the  different 
values  of  Z's  as  a  function  of  y  (x  held  constant,  x  >  f(^8) 
and  constructs  a  new  polynomial  z  «■  z  (y) .  For  the  same 
reason  as  before,  the  Lagrangian  form  was  used.  By 
interpolation  the  approximate  value  of  z  at  the  selected  y  - 
Yg  is  found.  This  is  the  value  of  Zg  «  Z(Xg,Yg)  at  the 
selected  X's  and  Y's. 

C.  THE  COMPUTER  CODE 
1.  gfiDggftl 

The  computer  code  was  written  originally  in  ST- 
BASIC.  The  code  includes  the  following  modules: 
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a)  START-’^startB  tha  program,  daflnoa  valuaa  of 

paramatara  and  gats  input 

b)  CDCURVB<— ‘daflnaa  tha  function  of  •  C^(efl^) 

c)  DENCRN— draws  tha  curvas  on  tha  scraan 

d)  CDCALC— •calcuiatas  tha  valua  of  by  doubla 

intarpolation 

a)  NawCon — prasants  tha  rasults  and  draws  tha  shape  of 
tha  new  cona 

f)  Weight— calculates  and  presents  astinatas  of  the  cone 
weight  and  volume. 

2.  Flow  Charts 

Flow  charts  of  the  program  code  are  shown  in  Figure 

2.4. 


D.  USING  THE  MODEL 

1.  Input 

As  input  to  the  program  the  following  variables  are 

needed : 

a.  The  cone  half  angle 

b.  The  ratio  of  cona  nose  radius  to  base  radius. 

To  assist  the  user  the  input  screen  represents  the 
data  bas3  and  the  meaning  of  the  input,  graphically. 

2.  Output 

a.  General 

The  output  of  this  model  is  given  by  graphic 
representations  and  numerical  values  on  the  screen.  A 
typical  output  is  shown  i.i  Figure  2.5. 
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Figure  2.4  Flowchart  of  CD  Cone  Program 
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Figure  2. 


5  Typical  Output 
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b.  Graphics 


Ths  graphical  raprasantation  includas  tha 
following  alaaants: 

1)  Tha  graph  of  a  •  s(x,y) 

2)  on  tha  graph  tha  valuaa  of  tha  salactad  variablaa 

Ya  ara  shown  basida  tha  valua  of  tha  dapandant 
variabla  (2,  -  Z(Xa,Y„)). 

3)  A  scala  schaao  of  tha  body  is  drawn.  Tha  achaaa 
providas  tha  dasignar  with  a  graphic  raprasantation. 
Tha  location  of  tha  cantar  of  prassura  (X^.p,)  is 
sarkad  claarly  on  tha  syaaatrie  axis  of  tha  hodf* 

c.  Nunarical  Rasults 

Tha  valuas  of  tha  following  paranatars  ara  shown 
digitally  on  tha  output  scraan: 

1)  Coafficiant  of  drag  (C(ji) 

2)  Tha  cona  naas  (Np) 

3)  Tha  cona  voluaa  (Vp) . 

Estiaatad  valuas  of  tha  projactila  voluaa  and 
mass  ara  calculatad  if  dasirad.  Tha  calculation  is  based 
upon  tha  aaasura  of  tha  projactila  basa  radius. 


I 

1 

1 


Ill 


A.  GENERAL 

1.  Bactairound 

El^ctronagntttlc  rail  launchar  concapta  (rallguna) 
hava  axistad  slnoa  tha  aarly  1900*8.  Tha  work  of  Rashlaigh 
and  Marshall  (1978)  and  Barbar  (1972)  gava  cradanca  to  tha 
potantlal  of  rallgun  accalarators .  Tha  raaaarch  was  ravivad 
at  tha  Australian  National  Univarsity  (1972)  using  a  single 
larga  Inductor  as  tha  po%rar  source.  More  recently  a 
Ijawronca  Livermore  National  Laboratory/Los  Alamos  National 
Laboratory  team  has  successfully  demonstrated  tha  use  of 
explosively  driven  magnetic  flux  compression  generators  to 
power  a  variety  of  rail  launchers.  Theoretical  and 
practical  researches  and  works  are  continuing  in  several 
directions,  predominantly  in  the  United  States.  This  model 
is  based  mainly  on  tha  open  literature  about  railgun 
studies.  The  studies  sponsored  by  NASA  concerning  the 
feasibility  assessment  of  earth-to-space  railgun  launchers 
for  launching  nuclear  waste  disposal  into  space  were  found 
to  be  valuable  references.  They  were  also  used  to  examine 
and  validate  this  model. 

2.  Railgun— The  Basic  Principles 

An  electromagnetic  railgun  consists  of  two 
conducting  rails  (electrodes) ,  between  which  a  conducting 
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elemen'b  with  an  attached  projectile  is  placed.  Electric 
current  is  passed  along  one  electrode  through  the  conducting 
portion  (armature)  of  the.  conductor/projectile  and  back 
along  the  other  electrode  (see  Figure  3.1).  The  current,  I, 
flowing  through  the  projectile  armature,  interacts  with  the 
magnetic  flux  generated  by  the  current  loop,  resulting  in  an 
I  X  B  force  in  the  direction  parallel  to  the  rails  (as 
indicated  on  Figure  3.1).  The  projectile  is  free  to  slide 
along  the  rails.  The  X  x  B  force  accelerates  the  projectile 
as  long  as  the  current  continues  to  flow  and  the  conductor 
remains  in  electrical  contact  with  the  rails. 

3.  Power  Supply  and  Conditioning 

In  order  to  fire  the  projectile  from  the  railgun, 
power  has  to  be  supplied  from  a  power  source.  A 
conditioning  system  has  to  transfer  this  power  to  suitable 
wave  form  that  will  assure  maximum  efficiency. 

The  dema.id  from  the  power  supply  system  is  dependent 
both  on  the  total  input  energy  needed  per  shot  and  the  rate 
of  shooting.  The  power  supply  system  must  be  suitable  to 
the  carrier  of  a  railgun.  Basically  the  following  main 
options  for  power  supply  exist: 

a.  Static  ground  system — in  this  case  power  may  be 
supplied  either  from  the  local  power  grid  or  from  a 
special  system  designed  for  this  purpose. 

b.  Movable  terrain  systam — this  system  would  rely  on  a 
self-contained  power  system  that  would  assure 
continuous  operation  of  the  weapon  system. 

c.  Seaborne  system — in  this  case  the  power  supply  would 
originate  either  from  the  vessel  (ship/ submarine) 
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Figure  3.1  The  Railgun 
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power  system  or  from  a  special  power  supply  system 
dedicated  to  this  mission. 

d.  Airborne  system— -same  options  as  in  the  seaborne  case. 
In  addition^'  the  power  system  has  to  be  light  and 
comply  with  the  standards  and  demands  for  airborne 
systems  4 

e.  Space  system— power  supply  is  difficult  to  achieve  and 
expensive  to  produce.  The  magnitude  of  power  needed 
to  fire  a  railgun  would  require  an  enormous  solar 
panel  or  a  large  nuclear  power  source. 

The  power  supply  from  the  primary  source  is  not  in  a 
suitable  form  to  be  used  by  the  railgun.  A  power 
conditioning  system  is  needed  in  order  to  produce  the  huge 
pulse  of  power  for  such  a  short  duration. 

The  power  conditioning  system  is  demonstrated  in 
Figure  3.2.  The  power  from  the  primary  supply  unit  is 
transformed  by  a  power  transmission  bus  to  a  DC-UC  converter 
that  regulates  the  voltage  in  order  to  transfer  it  to  the 
energy  storage  device.  The  energy  is  kept  for  a  short 
period  in  the  energy  storage  unit  until  the  moment  when  the 
gun  is  fired  and  the  switch  is  closed. 

A  pulse  forming  network  creates  the  final  shape  of 
the  power  and  sends  it  to  I'he  rails. 

The  short-time  storage  contains  the  high  energy  just 
before  sending  it  in  its  suitable  form  to  the  railgun.  It 
may  be  either  an  electric  system  (capacitive  system  based  on 
energy  storage  by  capacitors)  or  an  inductive  system  (in 
which  the  energy  is  stored  in  a  mechanical  form  as  a 
rotational  speed  of  a  iieavy  flywheel  of  a  homopolar 
generator) . 
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ERROOUCCO  AT  GOV 


Figure  3,2  Power  Conditioning  Systeri 
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The  Inductive  based  systeies  are  much  more  efficient 
for  s'.oring  energy  and  producing  the  suitable  waveform  fbr 
the  rallgun.  Most  efforts  in  recent  years  were  aimed  at 
improving  those  systems  and  to  make  them  suitable  for 
airborne  and  space  operation.  The  main  interest  is  now 
focused  on  the  distributed  energy  railguns. 

In  the  distributed  energy'  store  railguns  the  energy 
sources  required  to  power  the  gun  are  distributed  along  the 
length  of  the  barrel.  When  firing,  a  traveling  wave  of 
current  accelerates  the  projectile  by  sequential  switching 
of  energy  sources  distributed  along  the  rails. 

The  first  major  demonstration  of  this  idea  was  the 
construction  of  the  homopoiar  generator  (HPG)  at  the 
Australian  National  University  (ANU)  in  Canberra.  This 
machine  demonstrated  that  it  was  possible  to  make  very  large 
electromechanical  energy  stores  in  the  range  of  1  Gj  range. 
The  energy  is  stored  in  the  form  of  rotational  kinetic 
energy  of  heavy  (ten  tons)  rotors.  When  needed  the  full 
energy  can  be  extracted  electrically  into  a  suitable 
circuit. 

It  was  also  demonstrated  by  the  Canberra  HPG  that 
solid  brushes  can  be  used  to  carry  the  huge  currents  for  the 
short  acceleration-time  involved. 

The  success  of  the  Canberra  HPG  led  other  groups  to 
apply  the  same  techniques. 
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At  the  University  of  Texas  at  Austin,  a  10  Mj 
machine  was  operated  in  1981.  It  was  found  that  these 

machines  were  not  efficient  enough  and  only  about  6  percent 
of  the  mass  of  the  energy  was  actually  used  for  energy 
storage.  This  situation  led  to  a  new  concept  in  inductive- 
based  power  conditioning  systems  in  which  all  the  mass  of 
the  magnetic  circuit  can  be  used  as  a  rotational  energy 
device.  This  model  is  based  on  the  most  recent  data  about 
these  machines  that  were  available  in  the  open  literature. 
Other  new  concepts  of  inductive  systems,  like  the  reverse 
railgun  concept  (fast  moving  armature  is  producing  the 
current  wave)  and  the  explosive  driven  flux  compression 
concept  were  not  modeled  in  this  work. 

B.  THE  MATHEMATICAL  MODEL 
!•  List  of  Constants 

The  list  of  constants  used  in  this  model  is  given  in 
Table  3.1.  The  following  description  is  given  for  each 
constant : 

a)  the  program  code  symbol 

b)  the  text  symbol 

c)  the  definition 

d)  the  value 

e)  the  units. 

2.  List  of  Variables 

The  list  of  variables  used  in  this  model  is  given  in 
Table  3.2.  The  following  is  given  for  each  variable: 
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TABLE  3.1 


LIST  OF  CONSTANTS 


SYMBOLS 

MEANING 

VALUES 

UNITS 

Computer 

Code 

Text 

PI 

TT 

Pi 

3.14 

Non  Dlmen 

F 

f 

Screen  ratio 

79/33 

Non  Dimen 

Lo 

Inductance 

3 

Uo 

Mag.  permeability 

1.26* 

10“® 

H 

Ur 

Mr 

Resistivity  coper. 

20*10“9 

Ohm ' s/m 

AMAX 

^max 

Max  acceleration 
of  the  projectile 

100,000 

g's 

SPCWP 

Specific  weight 
of  the  projectile 

3000 

KG/m3 

SPCWR 

Pr 

Specific  weight 
of  the  rails 

9000 

KG/m3 

SPCWS 

Specific  weight 
of  the  structure 

2000 

KG/m3 

RYIEIJ) 

Yield  stress  of 
the  rails 

600,000,000 

N/m2 

KGKWl 

%)A 

Specific  weight 
per  power — 
airborne 

0.2 

Kg/KWATT 

K6KW2 

Specific  weight 
per  power- 
maritime 

0.4 

Kg/KWATT 

KGKW3 

Specific  weight 
per  power— solar 

5.3 

Kg/KWATT 

EDC 

Efficiency  of  the 
DC-DC  converter 

0.92 

Non  Dlmen. 

EC 

Efficiency  of  the 
capacitive  storage 

0.97 

Non  Dimen. 
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TABLE  3.1  (CONTINUED) 


SYMBOLS 

MEANING 

VALUES 

UNITS 

Computer 

Code 

Text 

E 

^in 

Efficiency  of  the 
inductive  storage 

0.99 

Non  Dimen 

EL 

^L 

Efficiency  of  the 
pulse  forming 
network 

0.99 

Non  Dlmen 

Ehp 

^hp 

Efficiency  of  the 

honopolar 

generator 

0.8 

Non  Dimen 

TLN 

Tln 

Temperature  of 
liqpjid  nitrogen 

77 

©K 

TACM 

Max  rail  temp. 

45® 

Or 

Ps.l. 

Density  of  air  at 
sea  level 

1.25 

Kg/m3 

FAIR 

Fcl 

Ratio  of  cooling 
area  to  surface 
area 

10 

Non  Dimen 

TABLE  3.2 
LIST  OF  VARIABLES 


SYMBOLS  MEANING  UNITS 


Computer 

Cod' 

Text 

T 

t 

Time  from  the  start  of 
acceleration 

Sec 

X 

X 

X  coordinate  along  the 
barrel  length 

m 

VP 

V 

p 

Velocity  of  the  projectile 
(along  X) 

M/sec 

33 


TABLE  3.2  (CONTINUED) 


SYMBOLS 

MEANING 

UNITS 

Conputsr 

Codtt 

Taxt 

Vmax 

^nax 

Max  valocity  of  tha 
projactila  (along  X) 

L 

L 

Inductanca  par  unit 
langth  of  rails 

Henry 

Lp 

L* 

Inductanca  par  unit 
langth  of  rails 

Harry/n 

Lo 

Inductanca  of  tha  storage 
and  wiring 

Henry 

R 

R 

Raaistanca 

Ohns 

Lmax 

%ax 

Max  allowabla  value  of 
barrel  langth 

n 

RP 

Rp 

Raaistanca  per  unit 
length  of  rails 

Ohm/meter 

Ro 

1 

Ro 

Resistance  of  tha  storage 
and  wiring 

Ohm 

i  I 

I 

Electric  current 

Amp 

lo 

lo 

Electric  currant  tine  T  »  o 

Amp 

F 

F 

Force 

Newton 

E 

E 

Energy 

Joule 

Eras 

Er 

Resistance  energy 

Joule 

Eloss 

El 

Energy  loss 

Joule  ^ 

1 

E indue 

Ei 

Inductance  energy 

Joule  1 

DPI 

Epi 

Energy  for  power  inputing 

Joule 

Ek 

Ek 

Kinetic  energy  of  the 
projectile 

Joule 

Es 

Es 

Energy  storage 

Joule 
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TABLE  3.2  (CONTINUED) 


SYMBOLS 

MEANING 

UNITS 

Computer 

Code 

Text 

EFRrg 

e 

Efficiency  of  the  railgun 

Non  Dimen. 

EFc 

fe 

Efficiency  of  the  capac. 
energy  store 

Non  Dimen. 

EFpc 

epc 

Total  power  conditioning 
efficiency 

Non  Dimen. 

B 

B 

Coefficient  of  rail  drag 

K/(m/S) 

H 

h 

Height  of  rails 

m 

W 

a 

Distance  between  rails 

m 

Ar 

AR 

Rails  cross  area 

m^ 

dS 

ds 

The  structure  width 

m 

AS 

AS 

Structure  cross  area 

10^ 

MB 

Mb 

Mass  of  the  barrel 

Kg 

Mr 

Mj. 

Mass  of  the  rails 

Kg 

Npwrs 

Mpg 

Mass  of  the  power  supply 
system 

Kg 

Mpwrc 

Mpe 

Mass  of  the  total 
conditioning  system 

Kg 

Epc 

Epc 

Efficiency  of  conditioning 
system 

Non  Dimen. 

Mbc 

Mbc 

Mass  of  the  power 
transmission  bus 

Kg 

Me 

Me 

Mass  and  efficiency  of  the 
capac.  energy  stor. 

Kg 

Msw 

Msw 

Mass  and  efficiency  of  the 
switch  of  the  rails 

Kg 

Mdc 

Mdc 

Mass  and  efficiency  of  the 

Kg 

DC-DC  converter 
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TABLE  3.2  (CONTINOBD) 


SYMBOLS 

MEANING 

UNITS 

CoBiput«r 

Cod* 

Taxt 

Npf 

Ml 

Mass  and  affioiancy  of  tha 
pulaa  foraing  ayataa 

Kg 

Nhpg 

Maaa  of  tha  hoaiopolar 
ganarator 

Kg 

Mind 

Min 

Maaa  of  tha  induotanca 
povar  atoraga 

Kg 

MSH 

Msw 

Maaa  of  tha  awitching 
natvorlc 

Kg 

Mgnt 

Mg 

Maaa  of  tha  quantity  of 
anaunition 

Kg 

Mstro 

M.t 

Maaa  of  railgun  atructura 

Kg 

Mtotal 

Mt 

Total  naaa  of  tha  full 
railgun  ayataa 

Kg 

V 

V 

Voluaa 

a^ 

Vind 

Vin 

VolvLna  of  tha  induotanca 
powar  atoraga 

Vhpg 

Vhg 

Voluaa  of  tha  homopolar 
ganarator 

Vind 

^ind 

Voluaa  of  tha  induotanca 
powar  atoraga 

a^ 

VSW 

^sw 

Voluaa  of  tha  switching 
natwork 

a^ 

PPS 

f 

Kata  of  fire 

Pulse/sec 

PWR 

^wr 

Average  powar  supply  to 
rails 

Watts 

K6PKW 

Wp 

Specific  weight  per  power 

Kg/KW 

FE 

Fm 

Electroaagnetic  force 

N 

FF 

Ff 

Friction  force 

N 
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TABLE  3.2  (COHTIMUSD) 


SYMBOLS 

NBAMZNG 

UNITS 

3aput«r 

Cod« 

Taxt 

Mp 

Mp 

Total  Baas  (projaotila— ‘ 
sabot) 

Kg 

BSD 

£ 

Blaotrical  skin  dspth 

B 

ta 

Acoalsration  tisa 

sac 

R1 

Ri 

Long  radius  of  all ipsa 

B 

R2 

Ra 

Short  radius  of  all ipsa 

B 

TIMDIJC 

Tin 

Tamparatura  of  inductor 

OK 

TRAD 

Tr 

Taatparatura  of  radiator 

Ok 

TAZR 

Tair 

Tonparatura  of  air 

Ok 

TWTR 

Taapara^vixa  of  soa  watar 

Ok 

TR 

Tr 

Tanparatura  of  rails 

Ok 

DETEMP 

DT 

Tauparatura  diffaranca 

Ok 

FAIR 

Ph/Ps.l 

Danaity  ratio 

VAIR 

Vair 

Airspaad 

m/a 

MTOTAL 

Kt 

Coafficiant  of  haat 
(TOTAL) 

W/KOcm2 

KWTR 

^wtr 

Coafficiant  of  heat  transfer 
to  water 

W/KOcm2 

KCOPPER 

Kcp 

Coefficient  of  heat  of  copper 

W/KOcm2 

KAIRW 

Kalr 

Coefficient  of  heat  transfer 
to  air 

W/KOcm2 

KAIUM 

Kal 

Coefficient  of  heat  of 
aluminuB 

W/KOcm2 

TALUM 

Kai 

Temperature  of  aluminum  plates  °K 

TCOPPER 

Tcp 

Temperature  of  copper  plates 

Ok 
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TABLE  3.2  (COMTZNUSD) 


SYMBOLS  MEANING  UNITS 


Computer 

Text 

Coda 

NCOOL 

*«cl 

Maas  of  tha  cooling  system  Kg 

VCOLL 

Vcl 

Volume  of  tha  cooling  system  m^ 

a)  Th«  prograa  coda  syBbol 

b)  Th«  taxt  sysibol 

c)  Tha  dafinition 

d)  Tha  units. 

3.  Basic  Aanunptlons  and  Ralationa 
a.  Ganaral 

Tha  assumptions  and  ralationa  that  ara  uaad  to 
conduct  this  modal  ara  dascribad  balov.  Thay  vara  gatharad 
in  groups.  Each  group  is  ralatad  to  ona  of  tha  different 
physical  aspects  of  the  modal.  Tha  dascription  of  the 
following  groups  of  assximptions  is  given  in  this  chapter: 

1)  Tha  railgun  modal 

2)  The  railgun  structure  and  mass  estimation 

3)  Tha  power  supply 

4)  Capacitive  based  power  conditioning  system 

5)  Inductive  based  power  conditioning  system 

6)  The  cooling  system 

7)  The  system  efficiency. 
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b.  Th«  Rftilgun  Nod«l 

(1)  Th#  Kauation  of  Motion.  Th«  •l«otroaag- 
n«tlo  forea  (Fsm)  acting  on  tha  projaotiia  la  givan  by  [Raf. 

41s 

Fbi  -  L'lVa 

vhara: 

L*  ia  tha  inductanoa  par  unit  langth  of  tha  rails 
I  is  tha  currant. 

Tha  aquation  of  notion  of  tha  projactila  in 

tha  barral  is: 


Frh  -  ^f 


vhara: 

N  is  tha  mass  of  tha  projact  (plus  tha  sabot) 
is  tha  projactila  valocity 
Fqi  is  tha  alactroaagnatic  forca 
Ff  is  tha  friction  and  drag  forca. 

Tha  alactroaagnatic  forca  in  tha  railgun  is 
rising  to  huga  valuas.  In  ordar  to  achiava  good 

conductivity  with  low  friction  tha  projactila  is  aounted 
during  tha  accalaration  process  in  a  specially  designed 

sabot.  If  tha  railgun  is  located  in  low  altitude  (on  a  ship 

on  land  or  on  a  low  flying  airplane)  and  the  dense  air  may 
cause  aerodynamic  drag  on  the  projectile,  a  specially 
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designed  pumping  system  ma>  be  installed  in  order  to  achieve 
a  low  pressure  at  the  barrel  before  shooting.  Therefore,  it 
may  be  assumed  that  the  electromagnetic  accelerating  force 
is  much  larger  than  the  friction  force.  As  a  result,  the 
friction  and  drag  forces  may  be  considered  negligible  on 
this  model: 

^f 


Under  these  assumptions  the  equation  of 
motion  has  the  following  form: 


M 


L 


l2 


The  mass,  M,  is  the  total  mass  accelerated 
by  the  railgun.  It  includes  the  projectile  and  the  sabot. 
For  the  purpose  of  this  model  it  was  assumed  that  the  sabot 
mass  is  about  25%  of  the  projectile  mass.  [Ref.  5]  Once 
the  projectile  leaves  the  muzzle,  the  sabot  is  detached  from 
the  projectile  and  decelerates  because  of  its  high 
aerodynamic  drag. 

(2)  The  Electrical  Circuit.  The  electrical 
circuit  of  the  railgun  is  shown  schematically  in  Figure  3.3. 
This  circuit  may  be  presented  in  the  form: 

V  -  isr(LI)  -  RI  *  0 
dt 

The  Inductance  L  may  be  represented  as  a  function  of  the 
proji.  :tile  location  along  the  barrel  (x)  . 
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L  -  L'X 


By  using  this  relation  we  get: 

I 

\ 

Ski  +  L  1^  i 

dt  dt 

L'l  ^  +  L  ^ 
dt  dt 

L'I\|,+  L  dl/dt 
and  the  circuit  equation  now  has  the  form: 

V  -  L'lVp-  Lf  -  R1  .  0  . 

(3)  The  Railoun  Cycle.  In  order  to  fire  a 
series  of  projectiles  the  railgun  must  be  operated 
cyclically.  The  cycle  involves:  charging  of  intermediate 
storage  devices  by  the  main  power  supply,  loading  a 
projectile,  accelerating  it  through  the  barrel,  and  finally 
restoring  the  system  to  the  initial  conditions. 

This  firing  cycle  is  repeated  at  a  rate 
determined  by  the  rate  of  fire.  In  this  model  it  is  assumed 
that  an  analysis  of  a  single  cycle  may  represent  all  cycles 
in  a  series  of  repetitive  shots. 

The  process  of  transferring  the  electrical 
power  generated  by  the  primary  source  to  kinetic  energy  of 
the  projectile  is  involved  with  considerable  energy  losses. 
The  main  mechanisms  and  magnitudes  of  energy  losses  during 
the  process  were  modeled  here. 
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For  an  inductive  driven  railgun  the  process 
starts  by  accelerating  tho  homopclar  generator  by  an 
electric  motor,  a  gas  turbine  or  any  other  device.  During 
the  charging  of  the  inductor  energy  is  dissipated  in  the 
inductive  store  and  the  charging  circuits  due  to  resistive 
losses.  When  the  switch  is  operated  current  flows  into  the 
accelerator.  Most  of  the  energy  is  Imparted  into  the 
pellet,  the  rest  is  dissipated  resistively  in  the  rails. 
After  the  process  is  over,  some  residual  magnetic  energy  may 
still  be  stored  in  the  inductor  and  the  accelerator.  A 
large  fraction  of  this  energy  may  be  dissipated  by  arcing  at 
the  muzzle.  The  energy  is  deposited  in  the  thermal  and 
kinetic  energy  of  plasma. 

(4)  The  Resistive  Losses.  The  resistive  losses 
during  a  cycle  of  acceleration  may  be  divided  into  two 
parts:  circuit  resistance  and  rail  resistance.  The 
resistance  of  the  circuit  and  therefore  the  losses  can  be 
reduced  by  Increasing  the  cross  sectional  area  of  the 
conductor  or  by  use  of  the  superconductors  techniques.  For 
the  purpose  of  this  model  it  is  assumed  that  circuit 
resistive  losses  may  be  considered  negligible  compared  to 
the  rail  resistance  losses.  [Ref.  6] 

The  current  flows  in  the  rails  during  the 
acceleration  phase.  The  process  is  very  fast  and  the 
resistance  of  the  rails  changes  continuously  during 
acceleration  as  the  current  propagates  through  the  rails. 
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The  real stance  is  further  affected  by  the  "velocity  skin 
effect"  which  confines  the  current  to  a  thin  sheet  on  the 
rail  surface.  [Ref.  4]  The  finite  rate  of  current 
diffusion  into  the  rails  must  be  considered  in  order  to 
properly  evaluate  the  rail  resistance.  The  incremental 
resistance  equation  for  the  rails  is: 


2  nrVA<S 


where: 

Rj.  is  the  rail  resistance 
h  is  the  rail  resistivity. 

is  the  projectile  velocity 
h  is  the  height  of  the  rails 
(S  is  the  electrical  skin  depth. 


The  electrical  skin  depth  (6)  was  derived 
[Ref.  7]  from  the  diffusions  of  a  step  function  of  a  current 
into  a  conductor  as: 


5  ^ 

where : 

t  is  the  time  after  the  start  of  application  of  the 
current . 

The  equations  were  combined  to  give  the 
rail  resistance  as  a  function  of  time. 
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V 

Rr  -  (2(uourA)^^Vh)  /  ^  <Jt  . 

0 

(5)  Inductive  Losses.  The  Inductive  energy 
stored  in  the  circuit  is  given  by: 

Bind  -  +  L'x)i2/2 

where : 

Lq  is  the  storage  inductance 

L'  is  the  inductance  per  unit  length  of  the  accelerator 
X  is  the  position  of  the  projectile 
I  is  the  current. 

At  the  beginning  of  the  process  (X  =  0)  all 
the  energy  Is  stored  in  the  storage  inductor 

Eind(X-O)  “  ^  1V2 

During  the  acceleration  process  x  increases 
and  some  energy  begins  to  be  stored  inductively  in  the  rail. 
At  the  end  of  the  process  the  residual  energy  in  the  storage 
Inductor  is  conserved.  The  inductive  energy  on  the  rails  is 
lost.  As  a  result  the  loss  of  inductive  energy  per  cycle 
may  be  considered  to  be: 

Rind  =  L'xi2/2 

c.  Railgun  Mass  Estimation 

The  estimation  of  the  railgun  mass  is  based  on  a 
barrel  configuration  as  illustrated  in  Figure  3.4.  It  is  a 
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Figure  3.4  The  Barrel  Cross  Section 


1 


simple  rectangular  copper  bar  embedded  In  a  plastic 
structure  reinforced  by  fiber.  The  plastic  cover  serves 
also  as  an  Insulation  material  for  the  barrel. 

The  cross  section  area  of  the  rails  (Aj.)  Is 

given  by: 


Aj.  =  2h  dr 

Rail  thickness  (dr)  must  be  at  least  equal  to 
the  electric  skin  depth  of  copper  for  the  acceleration  time 
().  Thicker  rails  do  not  result  In  lower  resistive  losses! 
From  weight  consideration  In  this  model  the  rail  thickness 
was  assumed  to  be  the  minimum  necessary — equal  to  ohe  skin 
depth. 


I 
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Thft  mass  of  tha  ralla  depends  on  their  cross  | 

section  (Ax*)  r  their  density  (px>)  and  their  length  (X): 

Mr  -  Pr  Ax  X 

The  pressure  stress  on  the  barrel  is  carried 
mainly  by  the  barrel  structure.  From  conservative 
considerations  it  is  assumed  that  the  full  recoil  force  is 
acting  on  the  fiber  reinforced  plastic  in  which  the  rails 
are  embedded. 

The  typical  yield  stress  (a  g)  of  such  materials 
is  about  600,000,000  N/m^.  The  recoil  force  (Fr)  is  given 
by: 

Fr  -  M*a 

The  minimal  cross  section  area  (Ag)  needed  to 
withstand  the  load  is  given  by: 

As  *  ^r/^s  *  M*a/as 

The  mass  of  the  structure  (Mg)  is  given  by: 

Mg  =  Pg  Ag  X 

1 

j 

where  pg  is  the  density  of  the  fiber-reinforced  plastic 
(typical  value  of  density  might  be  around  3000  Kg/m^) . 

' 

The  mass  of  the  total  barrel  (M3)  is,  therefore, 
the  sum  of  the  "^ail^  lass  (Mr)  and  the  structure  mass  (Mg)  :  ! 

M3  —  Mr  Mg  . 
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A  v«ry  important  constraint  for  tha  design  of 
the  system  is  imposed  by  the  maximum  value  of  acceleration 
in  which  the  projectile  and  its  internal  systems  can 
withstand  safely.  The  model  does  include  a  check  on  this 
value  wherein  the  calculated  value  of  the  acceleration 
exceeds  the  maximxim  value  permitted,  the  model  sends  a 
warning. 

The  hej  jht  of  the  rails  is  primarily  a  function 
of  the  dimensions  of  the  projectile.  For  this  purpose  it 
was  assumed  that  the  mass  of  the  projectile  plus  the  sabot 
may  be  represented  by  a  prism  with  a  square  cross  section. 
The  length  of  the  prism  is  about  five  times  the  side  of  the 
square  section. 

If  the  average  density  of  the  accelerated  mass 
is  (typical  value  of  1500  Kg/m^)  is  used,  the  total  mass 
is  given  by: 

^a  *  ^p  ^s  “  Pa  ^a  ®  Pa  ®  ^r^ 

where : 

is  the  total  mass  accelerated 
Hp  is  the  projectile  mass 
Mg  is  the  sabot  mass 
hj.  is  the  rail  height. 

The  value  of  the  rail  height  is  therefore  given 
by: 
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hr  •*  [{Mp  +  M.)/5  ajVS 


If  an  alllptlcal  croaa  saction  is  chosan  for  tha 
barral,  tha  nat  structural  araa  of  tha  barrol  cross  saction 
is  givan  by: 


A0  »  TrR3^R2  -  hr  (2dr  +  w) 


whara : 

radii  of  tha  all ipsa 
hr  is  tha  rail  haight 

dr  is  tha  rail  width 

V  is  the  distance  between  rails. 


d.  Power  Supply 

The  modal  of  the  power  supply  system  is  based 

upon  the  following  assumptions  and  relations; 

(1)  The  total  power  supply  needed  is  a  function  of  Joth 
the  amount  of  energy  needed  per  cycle  and  the  number 
of  cycles  per  time  unit  (rate  of  fire) : 


PWR  -  Ejc  *  f/  E 


where : 

PWR  is  tha  energy  per  cycle 
f  is  the  rate  of  fire  (pulses  per  second) 
e  is  the  overall  efficiency  of  the  process. 

(2)  The  weight  of  the  power  supply  system  is  dependent  on 
the  power  rec[uired  and  the  type  of  the  power  supply 
system.  Typical  values  of  specific  weight  per  power 
(Wp)  were  chosen  to  represent  the  following  types  of 
systems : 

(a)  For  airborne  systems  (f2)“0.5  Kg/KWatts 

(b)  For  seaborne  systems  (f^) — 1.0  Kg/KWatts 
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(c)  For  solar  spatial  basad  systsns  (fg)— 5.3 
Kg/KHatts 

Land  bassd  systems  may  get  their  supply  either  from 
military  type  generators  or  from  a  local  electric 
grid. 

The  weight  of  nuclear  spatial  based  power 
systems  is  a  more  complex  function  of  the  power  produced. 
The  specific  weight  per  power  of  the  system  is  given  by 
[Ref.  7]: 

Wp  -  5.3  -  (PWR/1000  -  2000)  *  17  [Kg/KW] 

The  value  of  the  specific  weight  per  power  does  not  decrease 
beyond  the  value  of  1.0  [Kg/iCw]  no  matter  how  high  the  value 
of  PWR  will  be. 

e.  Capacitive  Based  Power  Conditioning  System 

The  mass  and  volume  of  the  capacitive  based 
power  conditioning  are  the  sum  of  masses  and  volumes  of  the 
major  elements.  The  models  of  the  subsystem  are  based  on 
[Ref.  8  ].  The  main  parameters  that  affect  the  size  of  the 
power  conditioning  system  are:  the  stored  energy  (per 
shot) ,  the  power  needed,  the  various  efficiencies  of  the 
subsystems,  and  the  characteristic  acceleration  time. 

The  capacitive  based  system  mass  is  therefore: 

Mpc  *•  Mgc  Mqc  +  Mq  +  Ml  +  Msw 

where : 

Mbc  =  7.7  *  10"®  pV2 

=  10  +  0.057  Ppc^/^  +  0.014  Ppc^/^  +  1  *  lO"^  Pp^ 
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Me  '  0.011  Bs 

Msw  ■  Eg/  g 

ML  ■  3.08  *  10“^  Eg 

Thtt  overall  capacitiva  baaad  ayatan  afflclancy 
la  glvan  by  tha  ralatlona  [Raf.  8]: 

^pc  “  ^  *  L 

whara: 

esc  ■  1  -  *66  *  10“®  pV2 
Cdc  **  0.082 
Cq  -  0,97 
■  0.99 

whara: 

Mpc#  cpc  tha  total  powar  conditioning  mass  and 

afficiancy 

Ngc*  powar  tranamiaaion  bua  mass  and  efficiency 

Noc»  epc  the  DC-DC  converter  mass  and  efficiency 

£C  the  capacitive  energy  store  mass  and 

efficiency 

Mx,f  the  pulse  forming  network  mass  and 

efficiency 

M5  switch  mass  and  efficiency 

P  is  the  total  average  power  output  of  the  source 

Pc  *  ^BC  P 

Pdc  is  the  input  power  to  the  DC  converter 

Eg  is  the  stored  energy  (for  one  shot) 


f.  Inductive  Bas«d  Poirar  SyatM 

Tha  syatan  concapt  is  llluatratad  In  Figura  3.5. 
A  gaa  turblna  waa  ehoaan  aa  tha  primary  powar  aourca  for 
thia  nodal.  Tha  turbina  shaft  ia  connactad  to  tha  honopolar 
ganarator  which  ohargaa  an  induotiva  anargy  atoraga.  A 
rapatitiva  awitch  ia  uaad  to  connact  and  diaconnact  tha 
railgun  from  tha  induotiva  anargy  atoraga  and  tha  ganarator 
circuit. 

Tha  turblna  night  ba  a  typical  turbo~ahaft 
angina  (lika  thoaa  uaad  in  turboprop  anginaa  or  hallcoptara) 
or  a  apacially  daaignad  turbina  activated  by  a  gas 
ganarator. 

Tha  nodal ing  of  tha  inductive  baaed  powar 
conditioning  syatan  ia  founded  nainly  on  recant  infomation 
gathered  from  publications  and  tha  literature.  Most  data 
came  fron  tha  IEEE  Transactions  on  Macmatics  dedicated  to 
tha  third  synposiun  on  Electronagnatic  launch  technology 
[Ref.  8]. 

In  thia  nodal,  the  volune  and  the  weight  of  the 
system  were  estimated  to  be  as  follows:  the  system  volume 
and  mass  are  the  sum  of  its  three  major  components: 

a)  the  homopolar  generators, 

b)  the  inductor,  and 

c)  the  switching  network. 
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Figure  3.5  Inductive  Power  System 
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Hence : 


M  ■  Mhp  +  Mjn  +  Mgw 
V  -  Vhp  +  Viji  +  Vgy 


where: 

^HP  mass  and  volume  of  the  homopolar 

generator 

MiNr  Vj)f  are  the  mass  and  volume  of  the  inductance 

Nsw»  mass  and  volume  of  the  switching 

network. 

The  volume  and  mass  of  the  homopolar  generation 
are  proportional  to  the  energy  stored  (for  one  shot). 

Mhp  *  Kjjhp  Es  “  Ejjhp  hp 
Vhp  =  Kvhp  Eg  =  Kvhp  E)^/ 

where : 

Eg  is  the  maximum  energy  storage  (per  shot) 

Emhp  specific  volume  (per  energy  unit) 

Kvhp  specific  volume  (per  energy  unit) 

£|^  is  the  efficiency  of  the  homopolar  generator. 

The  values  of  and  Kvhpg  fEsf*  8  3  are: 

Khhp  “  0.178  ton/mg  (10  ion/56  mg) 

Kyjip  =  0.068  cubic  ro/mg  (0.9^  hs/56  mg) 
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The  estimate  for  the  volxime  and  the  mass  of  the 
inductor  system  is  based  on  [Ref.  8].  The  estimate  is  done 
for  inductors  of  the  coaxial  type.  These  inductors  are 
chosen  because  they  produce  no  external  magnetic  fields  that 
may  damage  the  conveyor  (orbiter,  airplane,  ship  or 
submarine) .  Their  mass  and  volume  are  dependent  linearly  on 
their  temperature.  A  48  Mg  aliuninum  inductor  cooled  by 
liquid  nitrogen  has  a  mass  of  1.5  tons  and  its  dimensions 
are:  1.8  m — diameter  and  1.5  m — length.  At  151°  room 
temperature  (15°  C)  the  same  amount  of  energy  may  be  stored 
in  inductor  having  10  times  more  mass  and  about  twice  the 
size.  The  temperature  dependent  model  for  estimating  the 
inductor  mass  and  volume  is  therefore: 

"I"  “ 

ViN  =  [w(0.9)2i.5][1  +  7(Tin--TLN)) 


Where : 

Mju  is  the  inductor  mass  [Kg] 

VjN  is  the  inductor  volume  [cubic  meter] 

Tin  i®  't*'®  inductor  temperature  [Kelvin] 

Tln  is  the  liquid  nitrogen  boiling  temperature. 

The  bus  and  the  switching  network  of  the 
inductive  sysrem  are  about  one  to  two  orders  of  magnitude 
smaller  than  similar  devices  of  a  capacitive  based  [Ref.  8] 
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conditioning  system.  Therefore  a  good  estimate  for  their 
mass  and  volume  is  given  by: 

Msw  -  [1-74  *  10-fi  Ks/  'a3/30 

Vsv  ^  ^sv/^sv 

7,500  Kg/m 

where: 

Msv  is  mass  of  the  bus  and  the  switching  network 
is  the  volume  of  the  bus  and  the  switching  network 

Sgv  is  the  specific  weight  of  the  bus  and  the  switching 
network. 

g.  The  Cooling  System 

The  low  efficiency  of  the  railgun  system  implies 
that  large  amounts  of  wasted  heat  must  be  rejected.  For  an 
efficient  operation  and  in  order  to  prevent  structural 
problems,  the  rails  temperature  must  be  kept  under  450°  K 
[Ref.  8].  A  cooling  system  is  needed  in  order  to  reject  the 
wasted  heat  and  to  prevent  the  temperature  from  rising  above 
the  permitted  limit. 

For  space-based  systems  the  most  promising 
cooling  device  is  derived  from  the  heat  pipe  concept.  The 
estimated  mass  of  the  heat  pipe  system  is  the  following 
function  of  the  cooling  power  and  the  radiator  temperature: 

Mr  =  1.95  *  10®  Tr^^  Pj. 
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where: 


Mj.  is  the  mass  of  the  heat  pipe  system 
Tj.  is  the  temperature  of  the  radiator 
Pf.  is  the  radiated  power. 

To  estimate  the  size  and  the  weight  of  an 
airborne  cooling  system,  the  standard  methods  of  heat 
transfer  calculation  were  applied  [Ref.  8].  The  air 

approaches  the  cooling  fins  in  its  ambient  temperature  and 
density.  The  wasted  heat  flows  through  the  fins  and  is 
transferred  to  the  passing  air.  This  model  accounts  for  the 
following  relevant  parameters:  airspeed,  altitude,  fin  size 
and  fin  thermal  characteristics. 

The  overall  conductance  of  the  heat  (Kt)  has  the 
following  relation  with  the  metal  conductance  (Kal)  and  the 
metal-to-air  conductance  (Kaij.)  • 


The  value  of  the  metal-to-air  conductance  is  the 
following  function  of  the  airspeed  and  relative  density 
affected  by  the  altitude: 

Kair  “  [4  +  80/500(Vair-l00)](Pair/Ps.l.) 

The  ratlgun  temperature  is  limited  to  450°  K. 
The  ambient  temperature  of  the  air  is  a  function  of  the 
altitude.  The  same  model  of  the  atmosphere  as  for  the 
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flight  simulation  was  ussd.  The  difference  between  these 
two  temperatures  is  therefore: 


DT  -  Tr  -  Tair 


where: 

lyr  is  the  temperature  difference 

Tair  is  ambient  temperature  of  the  air 
Tair  *  T||ir(^) 

Trio  is  the  maximxua  rail  temperature  (450^  K) . 

The  estimate  of  the  surface  area  of  the  cooling 
system  is  given  by: 


Mci  -  *  Kci  *  Pci 


where : 

Mci  is  the  mass  of  the  cooling  surface 

Wci  is  the  average  width  of  the  cooling  surface 

pcX  is  the  specific  weight  of  the  fins  material 
(aliiminum  alloys) 

Aqi  is  the  cooling  system  "surface  area" 
where : 

K^ir  is  the  metal  to  air  conductance 
is  the  airspeed 

Pair/P s.c  i®  relative  density. 

The  amount  of  wasted  heat  is  a  function  of  the 
power  and  the  efficiency  of  the  whole  railgun  system. 
Therefore: 


Pci  • 


whttre: 

Pci  is  the  cooling  power 
P^j.  is  the  input  power 
e  is  the  railgun  efficiency. 

The  cooling  power  on  the  other  hand  is  given  by: 

Pci  -  Rr  *  DT  *  Aci  *  Fci 

where: 

Kip  is  the  heat  conductance 

DT  is  the  temperature  difference 

Aci  is  the  surface  area  of  the  cooling  system 

Fci  is  the  ratio  of  cooling  area  (including  the  fins) 
to  the  surface  area. 

The  estimate  of  the  size  and  the  weight  of  the 
seaborne  cooling  system  was  done  in  a  similar  way. 

The  following  differences  between  the  air- 
cooling  system  and  the  water  cooling  system  were  taken  into 
account: 

a)  The  heat  transfer  is  from  the  fins  to  the  flowing 
water  (dependent  on  the  ship  speed) 

b)  The  temperature  of  the  sea  (assumed  to  be  10°  C — 283° 
K) 

c)  The  fins  being  made  of  copper. 
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h.  Syst«m  Efficiency 

In  order  to  evaluate  the  overall  efficiency  ( e) 
of  the  accelerating  process,  an  efficiency  factor  is 
defined.  The  factor  is  the  ratio  of  the  kinetic  energy  of 
the  projectile  at  the  nuzzle  (E)()  to  the  input  energy  (Epj^) 
that  was  supplied  to  the  railgun  by  the  primary  power 
source .  Therefore : 


.  _  ^ 

The  input  energy  that  is  supplied  tc  the  system 
during  one  cycle  either  diverted  to  kinetic  energy  of  the 
projectile  or  is  lost  in  the  process.  So: 

®pi  “  Ek  ®1 

where : 

Epj^  is  the  total  energy  input 

E](  is  the  kinetic  energy  of  the  projectile  at  muzzle 

E^  is  the  energy  loss. 

4.  Method  of  Solution 

The  numerical  solution  of  the  differential  equations 
is  achieved  by  a  modified  Euler  method.  According  to  this 
method,  in  order  to  improve  the  accuracy  of  the  results,  the 
values  of  the  variables  are  taken  as  the  sum  of  their  value 
at  the  start  of  the  interval  plus  half  the  value  of  their 
differential  in  the  interval  dt.  According  to  this  method 
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thm  diff«r«ntial  equations  of  the  model  take  the  following 
form: 


AV  «  Li2  At/2  Mp 

AX  »  (V  +  AV/2)  t 

ARr  »  2(nrno/Tr(t+At/2))  V2(v+Av/2)At/h 

AI  -  -(L'(V+AV/2)  +  Ro  +  (Rn-*^Rr/2)  )lAt/(Lo+Li(X+AX/2) 

V  «  V  +  AV 

X  »  X  +  AX 

Rj»  *  ®r  ^  ^  ^r 

I  *  I  +  AI 

t  «  t  +  At. 

C.  THE  COMPUTER  CODE 

1.  (?.sngcal 

The  computer  code  was  written  in  ST-BASIC.  The  code 
includes  the  following  modules: 

a)  Input — defines  the  values  of  the  parameters  and 
receives  input 

b)  Simulation — numerical  solutions  of  the  model.  The 
simulation  stops  when  one  of  these  conditions  is  met, 
for  example  desired  speed  reached,  length  of  barrel 
exceeds  maxlmiun  value 

c)  Railgun — designs  and  presents  the  railgun  structure 
according  to  the  simulation  results 


61 


d)  Pow«r  syatMt— calculatM  and  praaants  the  main  ] 
paranatars  of  tha  alactrical  componants  (primary  power  I 
supply,  conditioning  system,  switches,  etc.) 

a)  Weight  and  volume-- gives  an  estimate  of  tha  weight  and 
volume  of  all  major  components 

f)  Cooling  system— gives  an  estimate  of  the  weight  and 

size  of  the  cooling  system.  The  estimate  is  based  on 
the  wasted  heat  power  and  on  the  environmental 
conditions  | 

j 

g)  Weight  and  volume— summarizes  the  estimates  of  the 
weight  and  size  of  all  major  components  of  the  system. 

The  sum  of  the  weight  and  the  volume  of  the  whole 
system  is  sho%m  numerically  on  the  screen. 

h)  End— this  module  controls  the  communication  with  the 
user.  Decisions  can  be  made  either  to  run  the  code 
again  or  to  chain  directly  to  the  flight  simulation 
program. 

2.  Flowcharts 

Flowcharts  of  the  program  code  are  given  in  Figure 


3.6. 


0.  USING  THE  MODEL 
1.  IDSI^ 

The  input  to  the  computer  includes  the  following 
parameters  and  initial  conditions: 

a)  The  projectile  mass  (HP) 

b)  The  initial  current  (Iq)  at  time  t  »  0 

c)  The  desired  velocity 

d)  The  maximum  value  of  the  barrel's  length  permitted. 

The  code  Includes  many  values  of  parameters.  These 
values  can  be  altered  by  standard  editing  of  the  code.  In 
addition,  the  cooling  system  module  aslcs  for  input  about  the 
environmental  conditions. 
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START 


.S^^nJLATIO^ 

YES 

END 


Figure  3.6  Flowchart  of  the  Railgan  Program 


2.  Output 


The  output  of  the  progran  code  is  given  on  the 
monitor  screen.  The  following  layouts  are  used: 

a)  Simulation— ‘the  dynamic  development  of  the 

acceleration  process  is  shown  both  graphically  and 
numerically  on  the  screen  (Figure  3.7). 

b)  Railgun  structure — the  main  design  parameters  based  on 
the  simulation  results  are  shown  both  graphically  and 
numerically  on  the  screen  (Figure  3.8). 

c)  Power  system  design — main  results  of  the  power  system 
estimates  are  shotm  numerically  on  the  screen  (Figure 
3.9) . 

d)  Cooling  system — the  estimates  of  the  cooling  system 
size  and  weight  are  shown  numerically  on  the  screen 
(Figure  3.10). 

e)  Power  conditioning*— main  results  of  the  power 

conditioning  system  are  shown  numerically  on  the 
screen  (Figure  3.11). 

f)  Weight  and  site— summary  of  the  weight  and  size  of  the 
main  components  and  the  whole  system  is  given  on  the 
screen  ( Figure  3.12). 
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Figure  3.7  Output  1 — The  Railgun  Simulation 


6b 
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REPRODUCED  Ai  GOVERNMENT  EXPENSE 


Figure  3.9  Output  3— The  Power  System 
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UBfeMk  M  n  Mfi  m-l 


REPRODUCED  AT  GOVERNMEj-JT  EXPENSE 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


Figure  3.11  Power  Conditioning 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


TBTftLliEISHTS 

THE  BARREL  SVSIEH 

i.995TflH 

THE  (imUNlTION 

3.750  TOH 

THE  dmONITlON  STORAGE 

0.938  TON 

THE  AmUNITlON  LOADER 

1.262  TON 

THE  POUER  COKDITIOHING 

1921.428  TON 

— 

THE  POUER  SUPPLY 

1S.421  TON 

THE  STRUCTURE 

194.573  rON 

THE  TOTAL  ^EISHT  IS 

2148.386  TON 

IS  SIHSLATIOH  EY/Hl  ?  | 

• 

Figure  3.12  Total  Weight  and  Size 


IV.  THE  FLIgHT.  glMVtATIPg  ^SP.PEil 

A.  GENERAL 

The  flight  simulation  program  is  the  core  of  this 
research.  It  was  written  in  order  to  simulate  a  full 
ballistic  flight  of  a  projectile.  The  program  treats 
separately  the  three  major  phases  of  a  ballistic  flight, 
namely: 

a)  The  launch 

b)  The  orbital  flight 

c)  The  reentry. 

The  input  to  the  program  includes  the  parameters  cf  the 
projectile  and  the  initial  conditions  at  launch, 

The  output  of  the  program  is  a  time  history  of  the  m?in 
variables  describing  the  flight.  The  output  may  be  given 
either  graphically  on  the  monitor  screen  or  as  listing  cf 
numerical  values  (either  on  the  screen  or  by  the  printer) . 

The  output  includes  the  numerical  values  of  the 
kinematical  variables,  the  aerodynamic  variables  and  the 
thermal  variables. 

The  three  phases  of  flight  (launch,  orbit  and  reentry) 
can  be  run  separately  in  order  to  permit  study  cf  any  one  of 
the  phases  of  flight.  Ordinarily  the  program  automatically 
changes  from  one  phase  to  another  when  necessary. 
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B.  THE  MATHEMATICAL  MODEL 
List  of  Constanta 

The  list  of  constant,  ussd  in  this  modal  is  givnn  in 

Tablo  4.1.  Tha  following  description  is  given  for  each 
constant : 


Tha  program  coda 

symbol 

b) 

The  text  aynbol 

c) 

Tha  daf inition 

<a) 

Tha  value 

e) 

The  units. 

TABLS 

LIST  OP  C0NSTAN7,'S 


SYMBOLS 

MEANING 

VALUES 

UNITS 

Computer 

Code 

Text 

RE 

Radius  of  earth 

6371 

Em 

PI 

•ir 

Pi 

3.14 

Ncn  Dimen. 

hb 

H 

Scale  height 

9200  , 

Em 

g 

g 

Gravity 

acceleration 

9,81 

m/  sec^ 

dQdM 

Heat  cf  Ablation 
Carbon-Carbon 

28 

Kj/kg 

DENHO 

Ps.l 

Air  density  (S<I,) 

1*25 

Kg/m3 

TCKO 

T3.I 

Temperature  ( S . L) 

15 

oc 

PHO 

Ps.l 

Pressure  (S,L) 

1 

atm. 

MUSL 

V’3.1 

Kinematic 
viscosity  (S.L) 

0.000158 

ft^/sec 

2.  List  Qg  Yacifthlgg 


The  list 

of  variables  used  in  this  model  is  given  in 

Table 

4.2. 

The 

following  is  given  for  each 

variable: 

a) 

The 

program  code  symbol 

b) 

The 

text  symbol 

c) 

The 

definition 

d) 

The 

units. 

TABLE  4.2 

LIST  OF  VARIABLES 

SYMBOLS 

MEANING 

UNITS 

Computer 

Code 

Text 

t 

t 

Time  from  the  start  of 
flight 

Sec 

xr 

X 

Coordinate  along  X  axis 

Meter 

z 

z 

Coordinate  along  y  axis 

Meter 

hr 

Altitude 

Meter 

Vr 

V 

Velocity  of  projectile 

m/sec 

deg 

Y 

Flight  angle  above  the 
horizon 

Degrees 

dar 

Y 

Flight  angle  above  the 
horizon 

Radians 

mp 

Md 

Mass  of  projectile 

Kgm 

rp 

Rp 

Radius  of  the  projectile 

m 

cd 

Cd 

Drag  coefficient 

Non  Dimen. 

ap 

Ap 

Frontal  area  of  the 
projectile 

.a2 
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TABLE  4.2  (CONTINUED) 


SYMBOLS 

MEANING 

UNITS 

computar 

Code 

Text 

beta 

3 

Ballistic  factor 

N/m2 

vrx 

Vrx 

X  component  of  velocity 

m/sec 

vry 

Vry 

Y  Component  of  velocity 

m/sec 

arx 

*rx 

X  component  of  acceleration 

m/sec^ 

ary 

*ry 

Y  component  of  velocity 

m/sec^ 

ARAD 

9 

Polar  coordinate  in  the 
orbit  plane 

Radians 

R 

R 

Radial  coordinate  in  the 
orbit  plane 

m 

Elevation  angle 

Radians 

DRIDTI 

dR/dt 

First  time  derivative  of  R 

m/sec 

DR2DT2 

d^R/dt^ 

Second  time  derivative  of  R 

m/sec2 

DAIDTI 

dA/dt 

First  time  derivative  of  A 

Rad/sec 

DA20T2 

d2A/d2t 

Second  time  derivative  of  A 

Rad/sec' 

DIP 

Zg 

Z  coordinate  of  surface 
of  earth 

km 

dqdt 

• 

q 

Rate  of  heating  per  unit 
area 

Watt/cm 

qh 

q 

Total  heating  per  unit 
area 

J/cm^ 

eroscm 

dl 

Depth  of  erosion 

cm 

TCH 

ThC 

Temperature  of  the  air 

Degrees 

Celsius 

TKH 

ThK 

Temperature  of  the  air 

Degrees 

Kelvin 
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TABLE  4.2  (CONTINUED) 


SYMBOLS 

MEANING 

UNITS 

Computer 

Code 

Text 

TNOSE 

To 

Temperature  at  the 
projectile  nos'^ 

Degrees 

Kelvin 

SOS 

a 

Speed  of  sound 

m/Sec 

MACH 

m 

Mach  number 

Non  Dimen. 

PNOSE 

Po 

Pressure  at  the  projectile 
nose 

Atm 

PH 

Ph 

Pressure  of  the  air 

Atm 

DENH 

Ph 

Density  of  the  air 

Kgm/m^ 

MU 

Kinematic  viscosity  of 
the  air 

ft^/sec 

LENGTH 

L 

Typical  length 

ft 

RNRB 

Pn/Pb 

Ratio  of  radius 

Non-Dimen, 

0 

0 

Aerodynamic  drag 

Newton 

S 

S 

Reference  area 

m^ 

dm 

dm 

Mass  erosion 

Kg 

ht 

ht 

Total  enthalpy 

j/Kgm 

hw 

hw 

Carbon-carbon  enthalpy 

j/Kgm 

3.  Basis  .AsamBBtioJis  and  R?lati9ns 

a.  General 

The  model  is  divided  into  three  main  parts: 
launch,  orbit  and  reentry.  The  basic  assumptions  and 
relations  in  each  part  are  described  separately. 
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b.  Launch 

(1)  Equations  of  Motion.  Classical  equations 
of  notion  were  used  in  this  model.  It  was  assumed  that  the 
projectile  is  moving  in  the  same  vertical  plane.  The  axes 
were  chosen  as  follows: 

(a)  X  axis-tangent  to  the  surface  of  earth  at  the  launch 
site. 

(b)  Z  axis — in  the  direction  of  the  local  vertical  at  the 
Instant  of  launch  and  pointing  upward. 

(2)  The  Earth.  A  simple  model  of  the  earth  was 
selected  in  which  the  earth  is  assxuned  to  be  a  perfect 
sphere,  with  a  radius  of  6,370  km.  It  was  assvuned  also  that 
the  mass  of  the  earth  is  equally  distributed  and  the 
gravitational  force  is  constant  and  pointing  into  the  (-z) 
direction. 

(3)  Zhfi _ Model  of  the  Atmosphere .  A  simple 

model  of  the  standard  atmosphere  was  selected.  According  to 
this  model  it  is  assumed  that  the  temperature  (T^)  is 
decreasing  linearly  from  T  ■  15°  C  at  sea  level  to  T  »  -56° 
C  at  an  altitude  of  36,000  ft.  At  higher  altitudes  the 
temperature  is  constant  with  a  value  of  T  »  -56°  c. 

It  was  assumed  also  that  the  relationship 
between  the  air  density  (pjj)  and  the  altitude  (h)  is 
exponential  and  may  be  represented  by  the  following 
function; 


Ph  “  Ps.l  *  exp(h/H)  . 


wherft  Pa.l  density  of  air  at  saa  l«v«l  (Pg.i  *  l>25 

Kgn/m^)  and  tha  scala  haight  (H)  is  9.2  KA.  This  model  was 
found  suitable  to  be  used  for  launch  and  reentry 
calculations  and  was  also  chosen  by  NASA  technical  paper 
2614  [Ref.  10]  for  similar  purposes. 

(4)  The  Aerodynamic _ ttadftl.  The  aerodynamic 

model  is  responsible  for  producing  the  values  of  the 
aerodynamic  forces  and  the  pressures  on  the  projectile. 

It  was  assumed  that  the  stagnation  pressure 
and  temperature  may  be  calculated  by  the  isentropic 
relations  for  perfect  gas: 

,  Y/(Y-x) 

To  -  Th*(1  +  (Y-1)/2*m2) 

5  Y/(Y-1) 

Po  «  Ph*(1  +  (Y-1)/2*m2) 

using  the  average  value  of  y  *  1.4. 

In  most  cases  the  projectile  velocity  is 
above  the  speed  of  sound.  In  this  case  it  is  assumed  that 
the  stagnation  pressure  on  the  nose  tip  is  related  to  the 
stagnation  pressure  before  the  shock  by  the  normal  shock 
relation  for  perfect  gas: 


l+((Y-l)/2)in^ 


(Y-1) 

(Y+1) 


1/1-Y 


Considering  the  fact  that  the  main  region 
of  interest  is  in  the  hypersonic  domain,  the  drag  on  the 


projttctiltt  was  calculated  by  using  the  relation: 


D  -  1/2  Ph  V2  S  Cq 

The  value  of  Cq  is  given  either  as  the  product  of  the  drag 
coefficient  prograia  or  as  a  direct  input  by  the  user. 

(5)  Heating  Analysis.  It  was  assumed  that  the 
projectile  is  coated  by  a  layer  of  ablation  material.  For 
this  case  a  well-established'  approximation  is  given  by  [Ref. 
10] : 

g  ■  Cp^  V**  \p 


At  tlie  stagnation  point: 

2  2  1/2 

,  (1-K).72  -[(l-H).72o)^-0.52o^] 

“  0.26 


where : 


.  tnol.  w^ght  of  air  . 

'nol.  of  the  ablative  material^ 


V4 


(h. 


-V 


a,  C,  N,  N  are  numerical  coefficients. 

Here  p  is  the  effective  heat  of  ablation 
per  unit  mass.  Therefore  the  mass  of  ablation  material 
required  per  unit  area  is  given  by: 

dm  «  q/ £ 

where  q  is  the  total  accumulated  heat  load  per  unit  area: 
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q  -  /  q  dt 


Knowing  tho  density  of  th«  ablativ* 
natarial  (p^)  tha  dapth  of  tha  aroaion  night  ba  calculatad 
by  tha  ralation: 


d  >  dm/p^ 

Tha  valuas  of  tha  coafficiants  C,  N  and  N 
ara  dafinad  by  tha  charactaristics  of  tha  flow  ai.  \  tha 
radius  of  tha  projactila.  For  P,  V  and  9ivan  in  mks 
units  tha  valua  of  C  is: 

C  -  (1  -  hw/ht) 

and  q  will  ba  in  Watt/cn^. 

For  laninar  convactlon  tha  valuas  of  tha 
constants  N  and  M  ara:  N  »  0.5  and  N  <■  3.0. 

For  turbulent  convection  the  value  of  M  is 
dependent  on  the  range  of  speed: 

(  3.7  for  V  >  4  kn/s 
M  »  < 

(  3.37  for  V  <  4  km/s 

and  N  »  0.8. 

The  characteristics  of  the  flow  are 
determined  by  the  Reynolds  number.  Reynolds  number  (Rf))  is 
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a  function  of  tha  valocity  (V) ,  typical  body  dlnanslona  (1) 
and  tha  klnaaatic  vlacoaity  (v) . 

A  siapliflad  aodal  for  tha  function  of  tha 
klnaaatlc  viscosity  vs.  altituda  was  salactad.  According  to 
this  modal  tha  valua  of  tha  kinamatic  viscosity  is  givan  by: 


whara 

^s.l  kinamatic  vlacoaity  at  saa  laval  (0.000158 

ft^/aac) 

PH^.l  ratio  batvaan  tha  danslty  at  tha  flight 

altituda  to  tha  danslty  at  saa  laval  (standard 
atmosphara) . 


In  this  varalon  of  tha  modal  it  was  assumad 
that  tha  ablativa  matarial  coating  is  composad  of  tha  state- 
of-tha-art  carbon-carbon  coat.  Tha  haating  and  danslty 
pararaatars  vara  chosan  accordingly. 

(6)  Tha  Ballistic  Factor .  Tha  ballistic  factor 
of  tha  projactila  was  dafinad  and  calculated.  The 
ballistic  factor  maasuras  tha  ratio  batvaan  tha  aerodynamic 
forces  and  tha  Inartia  of  tha  projectile.  Its  magnitude  in 
pressure  units  la  equal  to  tha  valua  of  tha  dynamic  pressure 
needed  to  decelerate  the  projectile  at  a  rate  of  ig  (9.82 
m/sec^) .  So; 


Ig  - 


l/2pV^SC^ 


aerodynanac  forces 
nass  of  projectile 


m 


or: 


0  ■  1/2  pV^  •“ 

It  was  axpsctad  that  tha  ballistic  factor 
(3)  would  hava  tha  sajor  affact  on  tha  parfomanca  of  tha 
projaetila  both  during  atnospharic  launch  and  ataospharic 
raantry. 

In  this  nodal,  as  tha  nks  unit  systam  was 
salactad,  tha  units  of  tha  ballistic  factor  (  3)  are 
Nawtons/m^ . 

c.  Orbit 

(1)  Tha  Equations  of  Motion.  Tha  aquations  of 
notion  for  tha  two  body  problam  in  spharical  coordinates  ara 
as  follows: 

r  -  r(92  cos^  <p*  ^2)  «  - 

••  ••  •  ♦ 
r  9  cos  +  2  ydcos  (Ji  -  2  Ye<^'  sin  ({»  »  0 


where 

r  is  the  distance  froa  tha  earth's  canter 
9  is  tha  aziauth  angle 
is  the  elevation  angle 
u  is  the  gravitational  constant. 

By  selecting  the  orbital  plane  as  the  plane 
of  reference,  the  elevation  angle  becomes  zero. 
Therefore: 
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lIMIIMir*****  ^  M  l/WM  kAal 


<t)  »  c 


and: 


cos  (j)  »  1 
sin  (p  -  0 

And  the  equations  of  motion  become: 

•  •  • 

r-reZ  a  . 

••  ••  • 

r6  +  2  re  «  0 


After  transferring  this  system  of 
differential  equations  to  the  equivalent  first  order  system 
it  becomes:  ^ 

V 

r  =  Vr 

6  =  U 

« 

y  =  -  y/Y^  + 

•  • 
fi  *  -  2  r/r 


This  form  Is  known  as  the  polar  form  of  the 
r  ell  method.  The  polar  form  of  the  classical  Cowell 
method  was  applied  in  this  model  purposely.  This  selection 
permits  ue  use  of  larger  integration  step  size  for  the  same 
truncation  error. 

For  greater  accuracy  the  step  size  (dt)  can 
be  reduced.  To  further  the  model  accuracy,  a  new  method  was 
developed.  Tnis  method  is  explained  i;:  the  appendix. 
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For  this  specific  model  it  was  assumed  that 

the  projectile  is  moving  in  the  same  solar  plane.  The  polar 

coordinates  were  chosen  so  that: 

a)  The  origin  is  at  the  center  of  the  earth. 

b}  The  angle  is  measured  from  the  launch  site  as  the 
reference  point  of  direction  (positive  in  the 

direction  of  the  projectile) . 

(2)  The  Earth.  Again  the  same  model  of  the 
earth  was  selected.  According  to  this  model  the  earth  is 
assumed  to  be  a  perfect  sphere  with  a  radius  of  6,370  km. 
The  mass  of  the  earth  is  concentrated  in  its  center. 

(3)  The  Space  Environment.  For  this  simplified 
model  it  was  assumed  that  the  density  of  air  in  space  is 
negligible  and  therefore  no  drag  is  encountered.  As  a 
result  the  total  energy  of  the  projectile  (the  sum  of  the 
kinetic  and  the  potential  energy)  is  constant  during  the 
whole  phase  of  orbital  flight. 

d .  Reentry 

The  reentry  model  simulates  the  flight  of  the 
projectile  from  the  moment  it  crosses  the  h  =  100  km 
altitude  line  and  enters  the  atmosphere. 

Basically  the  equations  are  similar  to  the 
equations  of  the  launch  program.  The  initial  conditions  for 
the  reentry  program  are  the  final  results  of  the  orbital 
module. 
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4.  Method  of  Solution 


In  all  three  modules  of  the  flight  simulation  the 
Euler  method  was  used  in  order  to  solve  the  differential 
equations. 

In  order  to  assure  optimal  ratio  of  accuracy  to 
running  time,  the  time  step  was  tailored  to  the  flight 
conditions.  During  the  launch  and  the  reentry  phases 
shorter  time  intervals  are  used,  during  the  orbital  flight 
phase  longer  time  Intervals  are  used  because  drag  effects 
are  negligible. 

In  addition  a  special  algorithm  was  developed  for 
solving  the  orbital  module.  This  algorithm  assures  accurate 
results  for  cases  when  longer  orbital  flight  is  simulated. 
This  method  is  explained  in  detail  in  Reference  11. 

C.  THE  COMPUTER  CODE 
1.  General 

The  flight  simulation  program  served  also  as  a  core 
for  some  other  programs.  The  following  programs  are 
essentially  based  on  the  same  mathematical  models. 

a .  Launch 

This  program  simulates  only  an  atmospheric 
launch  of  the  projectile  (altitude  less  than  100  km) . 

b.  Orbit 

This  program  simulates  only  an  orbital  phase  of 
a  flight.  It  assumes  zero-drag  and  therefore  can  be  used 
only  while  altitude  is  above  100  km. 
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c.  Reentry 

This  program  simulates  only  the  reentry  phase  of 
a  flight.  It  does  take  into  account  the  atmospheric  drag. 
Therefore  this  program  is  most  suitable  to  simulate  i^he 
reentry  from  an  altitude  of  100  km  and  belov. 

d.  Full  Flight  Simulation  (Summary  Printing 
Version) 

This  version  of  the  flight  simulation  program 
includes  an  additional  module  which  is  responsible  for 
automatically  sending  a  summary  of  the  results  to  the 
printer  at  the  end  of  each  run. 

The  flowchart  of  this  program  computer  code  is  shown 
in  Figure  4.1. 

D.  USING  THE  PROGRAM 
1.  Input 

The  input  to  the  flight  simulation  program  includes 
the  following; 

1.  Parameters  of  the  projectile: 

a.  Mass  of  the  projectile 

b.  Coefficient  of  drag  of  the  projectile 

c.  Radius  of  the  projectile 


2.  Initial  conditions: 

a .  Launch  speed 

b.  Angle  of  launch  (above  horizon) 

c.  Altitude  of  the  launcher  at  the  instant  of  firing. 
When  the  program  is  rerun,  the  last  set  of  initial 

conditions  and  parameters  is  used.  If  the  user  wishes,  he 
may  enter  a  new  set  of  values  or  change  a  selected  few  , 
separately. 

When  chained  (or  called)  by  the  previous  programs, 
such  as  the  drag  estimation  code  or  the  railgun  simulation, 
the  results  of  those  programs  are  passed  into  the  flight 
simulation  model  as  input. 

2.  QutPVt 

a.  General 

The  graphic  output  of  the  program  is  composed  of 
the  following  three  basic  layouts: 

1 .  Launch  layout 

2 .  Orbit  layout 

3.  Reentry  layout. 

When  the  printer  version  of  the  program  is  run 
either  a  listing  of  selected  variables  or  a  summary  of  the 
run  results  is  printed. 

The  summary  of  the  results  includes  the  values 
of  selected  variables  at  the  instant  of  the  following 
events : 
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1 .  Launch 

2.  Entrance  to  orbit  (H  ■  100  km) 

3.  Start  of  reentry  (H  »  100  km) 

4.  Ground  impact  (sea  level). 

b .  Launch  Layout 

An  example  of  the  launch  layout  is  shown  in 
Figure  4.2.  The  launch  layout  includes  graphic  representa¬ 
tions  and  nximerical  results.  The  time  history  of  these 
values  may  be  examined  by  the  user  as  he  follows  the 
development  of  the  graphs  and  the  changes  of  numerical 
values  during  the  simulation. 

The  graphic  representation  includes  the 
following  curves: 

1)  Speed  vs.  altitude 

2)  Deceleration  vs.  altitude 

3)  Rate  of  heating  vs.  altitude. 

A  picture  of  the  plane  of  flight  is  given  on  the 
other  side  of  the  screen.  The  contour  of  the  earth  is  shown 
under  the  projectile  location. 

The  values  of  the  following  parameters  are  shown 
numerically  on  the  screen  and  are  updated  every  second: 


1) 

t 

[sec] — the  elapsed  time  from  launch 

2) 

V 

[km/sec] — speed 

3) 

H 

[km] — altitude 

4) 

R 

[Nm] — ground  distance  from  launch  site 

5) 

g' 

s — Deceleration 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


Figure  4.2  Launch  Layout 
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M — mach  nuinbar 


6) 

7)  Tc — anblenc  tainparatura  [®C] 

8)  a — spaed  of  sound  [m/sac] 

9)  Tf|»stagnatlon  tanparatura  at  the  nose  tip 

10)  OQDT  [ kwatt/cn^ ] ••rata  of  heating  (at  nose  tip) 

11)  QH  [kg/cm^] — accumulated  heat  (at  nose  tip) 

12)  dMp  [q/cm^] ••loss  of  mass  due  to  thermal  erosion  (at 
nose  tip) 

13)  Eros  (cm) — depth  of  erosion  (on  the  carbon-carbon 
coated  nose  tip) . 

In  addition  the  following  parameters  of  the 
projectile  are  shown  constantly  on  the  screen: 

1)  — coefficient  of  drag 

2)  Mp — mass  of  the  projectile 

3)  Rp — radius  of  the  projectile  base 

4)  B— the  ballistic  coefficient  [N/m^]. 

c «  Orbit  Layout 

(1)  General .  The  orbital  layout  (Figure  4.3) 
includes  graphical  representation  and  numerical  results. 
The  time  history  of  the  results  may  be  examined  by  the  user 
as  he  follows  the  development  of  the  curves  and  the  changes 
of  the  numerical  results  during  the  simulation. 

(2)  The  Graphics.  The  graphic  representation 
shows  the  plane  of  motion  of  the  projectile.  The  earth  and 
its  center  are  shown  in  the  middle  of  the  screen.  The 
movement  of  the  projectile  is  shown  dynamically  on  the  same 
scale.  The  picture  is  updated  about  every  second. 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


Desk  File  Run  Edit  Debug 


Figure  4 . 3  Orbit  Layout 


“>* - ttUMtifial - EMuUa.  Th«  num«ric»l 

r«s«lt«  of  th*  following  v«riabl««  ,«  shown  dynamically  on 
ths  scraan  (thair  valuaa  ara  updatad  about  avary  aacond) : 

1)  t  C®«c]— «lapa«d  tina  from  launch 

2)  V  [km/s]  'Spaad 

3)  H  CJca]— altituda  (abova  aaa  laval) 

4)  R  [NM] —ground  distanca  from  launch  site. 

d.  Reentry  Layout 

A  typical  look  of  the  reentry  layout  is  shown  in 
Figure  4.4.  Basically  the  reentry  screen  is  similar  to  the 
launch  screen. 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


1 


Disk  File  Run  Edit  Debug 


Figure  4.4  Reentry  Layout 
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V.  THE  CTOVNO  mCK.«QPEL 


A.  GENERAL 

The  ground  track  model  is  based  mainly  on  the  orbit 
module  of  flight  simulation.  In  addition  to  the  orbit 
calculation  the  rotational  movement  of  the  earth  was  taken 
into  account.  The  results  produced  by  this  model  are 
represented  as  the  trace  of  satellite  orbit  over  the  world 
map  as  a  function  of  time. 

Input  to  the  model  includes  the  initial  condition  and 
the  parameters  of  the  orbit.  Output  of  the  program  includes 
a  graphic  picture  of  the  track  and  the  scheme  of  the  body 
movement  on  the  orbital  plane.  In  addition,  major 
characteristic  variables  of  the  flight  are  shown  numerically 
on  the  screen  or  are  printed  by  the  printer. 

B.  THE  MATHEMATICAL  MODEL 

1.  List  of-Congt.anl;§ 

The  list  of  constants  used  in  this  model  is  given  in 
Table  5.1.  The  following  description  is  given  for  each 
constant: 


a) 

The 

program  code  symbol 

b) 

The 

text  symbol 

c) 

The 

definition 

d) 

The 

value 

e) 

The 

units. 
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TABLE  5.1 


list  op  constants 


SYMBOLS 

MEfJHHG 

VALUES 

UNITS 

Computer 

Code 

Text 

. 

RE 

Re 

Radius  of  earth 

6371 

Km 

PI 

tr 

Pi 

3.14 

Non  Oimen 

hb 

H 

Scale  Height 

9200 

Kn 

g 

g 

Gravity 

acceleration 

9.81 

m/sec^ 

dQdN 

dq/dM 

Heat  of  Ablation 
Carbon-Carbon 

28 

Mj/kg 

OENHO 

^s.I 

Air  density  (s.l) 

1.25 

Kg/m^ 

TCHO 

Ts.l 

Temperature  (s.l) 

15 

Oc 

PHO 

Ps.l 

Pressure  (s.l) 

1 

atm 

MUSL 

^s.l 

Kinematic 
viscosity  (s.l) 

0.000158 

ftVsec 

2. 

List.  . 

yariablfla 

The  list 

of  variables  used  in 

this  model 

is  given 

Table  5.2.  The  following  is  given  for  each  variable: 


a) 

The 

program  code  symbol 

b) 

The 

text  symbol 

c) 

The 

definition 

d) 

The 

units. 
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TABLE  5.2 


LIST  OF  VARIABLES 

SYMBOLS 

MEANING 

UNITS 

Computer 

Code 

Text 

t 

t 

Time  from  the  start  of 
flight 

Sec 

VL 

Vl 

Speed  of  earth  rotation 

rad/ sec 

L 

dL 

Net  eastward  movement 

rad 

NL0N6 

New  longitude 

Degrees 

0L0N6 

k) 

Old  longitude 

Degrees 

TL0N6 

True  longitude 

Degrees 

DT 

dt 

Time  interval 

Sec 

Xr 

X 

Coordinate  along  X  axis 

m 

z 

z 

Coordinate  along  y  axis 

m 

hr 

Altitude 

m 

vr 

V 

Velocity  of  projectile 

m/sec 

d 

Y 

Flight  angle  above  the 
horizon 

Degrees 

dar 

Y 

Flight  angle  above  the 
horizon 

radians 

mp 

Md 

Mass  of  projectile 

Kgm 

rp 

Rp 

Radius  of  the  projectile 

m 

cd 

Cd 

Drag  coefficient 

Non  Dimen. 

ap 

Ap 

Frontal  area  of  trhe 
projectile 

m^ 

beta 

3 

Ballistic  factor 

N/m2 

vrx 

Vrx 

X  component  of  velocity 

m/sec 

96 


TABLE  5.2  (CONTINUED) 


SYMBOLS 

MEANING 

UNITS 

Computer 

Code 

Text 

vry 

Vry 

Y  component  of  velocity 

m/sec 

arx 

*rx 

X  component  of  acceleration 

m/sec^ 

ary 

*ry 

Y  component  of  velocity 

m/sec^ 

AEAD 

6 

Polar  coordinate  In  the 
orbit  plane 

Radians 

R 

R 

Radial  coordinate  In  the 
orbit  plane 

m 

4) 

Elevation  angle 

Radians 

DRIDTI 

dR/dt 

First  time  derivative  of  R 

m/sec 

DR2DT2 

Second  time  derivative  of  R 

m/sec2 

DAIDTI 

dVdt 

First  time  derivative  of  A 

Rad/ sec 

DA20T2 

dVdt^ 

Second  time  derivative  of  A 

Rad/sec 

DIP 

Ze 

Z  coordinate  of  surface 
of  earth 

km 

dqdt 

dq/dt 

Rate  of  heating  per  unit 
area 

watt/ cm 

qh 

q 

Total  heating  per  unit 
area 

J/cm^ 

eroscm 

dl 

Depth  of  erosion 

cm 

TCH 

ThC 

Temperature  of  the  air 

Degrees 

Celsius 

TKH 

ThK 

Temperature  of  the  air 

Degrees 

Kelvin 
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3 .  Baalc  Agaumptiona  and  Ralatlons 


The  nodal  is  based  mainly  on  the  orbit  module  of  the 
flight  simulation  program.  The  solution  of  Kepler's 
equations  was  achieved  by  Cowell's  method.  For  longer 
orbital  flights  or  when  nreater  accuracy  is  needed,  the 
solution  of  the  equations  is  achieved  by  an  energy  approach 
that  was  developed  as  part  of  this  research.  The 
description  of  this  method  is  given  in  Reference  11. 

Host  aastimptions  and  relations  described  in  Chapter 
II.A.3.C  are  valid.  In  addition,  the  rotation  of  the  earth 
was  Introduced.  Because  of  the  earth's  rotation,  th2 
eastward  angular  movement  of  a  body  with  respect  to  the 
earth's  surface  is  the  algebraic  difference  between  the 
eastward  angular  movement  of  the  body  with  respect  to  the 
center  of  the  earth  and  the  eastward  angular  movement  of  the 
surface  of  the  earth  (on  the  same  latitude)  with  respect  to 
the  same  reference  point. 

The  eastward  speed  (Vj^)  of  the  movement  of  the 
surface  of  the  earth  is  given  by: 

Vl  “  2Tr/ (3600*24)  [rad/sec] 

In  a  period  of  dt,  the  body  moves  from  the  longitude 
of  Lq  (old  longitude)  to  longitude  of  Ln  (new  l  ongitude) . 
The  net  movement  of  the  body  track  with  respect  to  the 
surface  of  the  earth  (the  point  on  the  surface  under  the 
body  location)  is  given  by: 
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dL  -  (Ln  -  Lo)  “ 


whttrs: 

dL  ia  th«  n«t  aastward  movamant 
li)l  Is  tha  naw  longltuda 
is  tha  old  longltuda 
Vx,  is  tha  body  aastward  velocity 
dt  is  tha  time  interval. 

4.  Method  of_  Solution 

Tha  sane  raathods  of  solution  as  for  the  orbit  module 
of  tha  flight  simulation  are  used  hare  to  solve  Kepler's 
equations. 

C.  THE  COMPUTER  CODE 

1.  General 

The  code  of  this  model  is  based  mainly  on  the  code 
of  the  orbit  module  of  the  flight  simulation.  In  addition 
to  the  solutlo  of  the  two  body  problem,  the  code  includes  a 
module  that  produces  the  world  map  on  the  screen  and  a 
second  module  that  locates  the  projectile  trace  on  the  same 
map. 

2.  Flowchart 

The  flowchart  of  this  program  is  shown  in  Figure 

5.1. 

3.  Basic  Code 

The  list  of  the  Basic  Code  is  given  in  Reference  11. 
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0.  USSiMG  THE  MODEL 
1.  Input 

Input  to  tho  program  includas  tha  following 
paramatars  and  initial  conditions: 

a)  Tha  initial  location  of  tha  body  in  aarth  coordinatas 
(latituda,  longituda,  altituda  [km]) 

b)  Tha  body  apaad  [km/sac] 

c)  Tha  body  climb  angla  [dagraas] 

d)  Tha  flight  diraction  [dagraas]  (flight  diraction  at 
tha  aquator  is  agual  to  tha  orbit  inclination) . 

A  typical  output  of  tha  program  is  shown  in  Figure 
5.2.  Tha  picture  on  tha  scraen  has  three  parts:  numerical 
results,  the  world  map  and  the  orbital  plane. 

The  values  of  the  following  variables  are  shown 
continuously  on  the  screen.  Their  value  is  updated  about 
every  second: 

a)  Tine  [sec]  is  the  elapsed  time 

b)  Velocity  [km/s] 

c)  Altitude  [km] 

d)  Latitude  [degrees] 

e)  Longitude  [degrees] 

In  addition,  the  values  of  the  following  parameters 
are  shown  as  well: 

a)  Longitude  of  A.N.  [degrees] 

b)  Initial  altitude  [km] 

c)  Inclination. 
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Figure  5.2  Typical  Output  of  the  Groundtrack  Program 
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Th«  ground  triico  of  tho  novlng  body  ovor  tho  surfaco 
of  tho  oarth  is  sho«m  continuously  as  ths  simulation 
advance.  Bach  cycle  is  shown  in  a  different  color. 

A  snail  schene  of  the  orbital  plane  is  shown  on  the 
upper  part  of  the  screen.  This  layout  pemits  the  user  to 
examine  the  shape  of  the  orbit  clearly. 
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VI .  THE  GUIDAMCE  MODEL 


A.  GENERAL 

Th«  yl«ld  of  a  warbaad  of  any  waapon  la  aainly  a 
function  of  ita  aaaa  and  voluna.  In  many  oaaaa  tha  ability 
to  anlarya  tha  warhaad  ia  liaitad.  Thia  ia  aapacially  trua 
in  tha  caaa  of  a  rallgun  whara  tha  aisa  and  naaa  of  tha 
projactila  hava  to  ba  kapt  aa  low  aa  poaaibla. 

In  ordar  to  incraaaa  tha  kill  probability  of  tha  anall 
railgun  projactila,  anothar  option  than  incraaaing  ita  aiaa 
has  to  ba  axaninad.  Tha  altarnativa  ia  to  achiava  graatar 
accuracy. 

Tha  goal  of  incraaaing  tha  accuracy  of  any  firing  systaa 
is  not  an  aaay  teak  to  achiava.  Tha  trajactory  of  tha 
projactila  ia  affactad  by  many  factora.  Evan  a  alight 
inaccuracy  in  tha  launch  conditions,  in  tha  projactila  siza 
or  in  tha  pradiction  of  tha  waathar  conditions  may  result  in 
a  significant  miss  of  tha  target. 

Increasing  tha  accuracy  of  tha  launching  system  is  not 
cost  effective.  Evan  with  considarabla  efforts,  only 
limited  improvement  in  accuracy  would  ba  achieved  and  the 
destruction  of  the  target  would  not  be  assured. 

The  tendency  therefore  is  to  provide  the  warhead  with 
some  final  homing  capability.  Homing  may  be  achieved  by  a 
passive,  a  semi-active  or  an  active  system.  It  is  more 
likely  that  railgun  projectiles  will  use  a  passive  system  as 
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thttir  apac*  and  wal^ht  ara  vairy  linitad.  In  thia  modal  a 
paaalva  homing  ayatam  waa  almulatad.  A  proportional 

guidanca  ayatam  with  variabla  conntanta  of  proportion  waa 
aalactad. 

Thia  modal  may  ba  uaad,  for  axampla,  for  aimulation  of  a 
warhaad  homing  tovarda  a  targat  illuminatad  by  a  lasar 
daaignator.  The  illuminator  may  ba  locatad  in  spaca  aboard 
an  airplana  or  on  land. 

B.  THE  MATHEMATICAL  MODEL 

1.  Lint  of  Conatanta 

Tha  liat  of  conatants  usad  in  thia  modal  is  given  in 
Tabla  6.1.  Tha  following  dascription  is  given  for  aach 
constant : 

a)  Tha  program  coda  symbol 

b)  Tha  taxt  symbol 

c)  Tha  dafinition 

d)  The  valua 

e)  Tha  units. 

2.  Lia^  Of.  Yarittbitta 

The  list  of  variables  usad  in  this  model  is  given  in 
Table  6.2.  The  following  is  given  for  each  variable: 


a) 

The 

program  code  symbol 

b) 

Tha 

taxt  symbol 

c) 

The 

definition 

d) 

The 

units. 
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TABU  6.1 
U8T  OP  CONSTANTS 


SYMBOU 

NBANINO 

VALUS8 

UNITS 

Computer 

Code 

Text 

BE 

Re 

Radius  of  earth 

6371 

Km 

PI 

ff 

Pi 

3.14 

Non  Dimen. 

hb 

H 

Soale  Height 

9200 

Km 

g 

g 

Gravity 

acceleration 

9.81 

m/sec^ 

dQdM 

dq/dn 

Heat  of  Ablation 
Carbon-Carbon 

28 

MJ/kg 

DEN-flO 

P«.l 

Air  density  (s.l) 

1.25 

Kg/m3 

TC+10 

Te.i 

Temperature  (s.l) 

15 

©C 

PHO 

P«.l 

Pressure  (s.l) 

1 

atm 

MUSL 

^•.1 

Kinematic 
viscosity  (s.l) 

0.000158 

ftVaec 

TABU  6.2 
LIST  OF  VARIABUS 


SYMBOU  MEANING 

Computer  Text 
Code 


UNITS 


t 

t 

Time  from  the  start  of 
flight 

Sec 

xr 

X 

Coordinate  along  X  axis 

Mater 

X  targ 

X^ 

Target  location 

Meter 

dX  targ 

dX 

Miss  distance 

Meter 
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TABLE  6.2  (COMTZIIUBD) 


SYMBOLS 

NBANZMO 

UNITS 

Co)lpUt«r 

Cod* 

T«xt 

ZRMTNBN 

M 

Vartical  baarino  of  tar^at 
(Haw) 

Radians 

ZRMTOLD 

0 

Vartical  baaring  of  targat 
(Old) 

Radians 

XOB 

K 

Constant  of  proportional 
navigation 

Non.  Diman 

NP 

‘'P 

Manauvaring  aocalaration 

g‘s 

COP 

Cd 

Coafficiant  of  drag 

Non  Oifflsn. 

CLP 

Cl 

coafficiant  of  lift 

Non  Dinan. 

COO 

Coo 

Coafficiant  of  paraaita 
drag 

Non  Dinan. 

K 

X 

Coafficiant  of  inducad 
drag 

Non  Dinan. 

CLP  Stall 

ClS 

coafficiant  of  lift  at 

Stall 

Non  Dinan. 

z 

z 

Coordinata  along  y  axis 

Matar 

hr 

Altituda 

Mater 

vr 

V 

Velocity  of  projactila 

Vaac 

d 

Y 

Diva  angla 

Dagraaa 

dar 

Y 

Diva  angla 

Radians 

np 

Md 

Naas  of  projactila 

Kgn 

rp 

Rp 

Radiua  of  tha  projactila 

n 

ap 

Ap 

Frontal  araa  of  tha 
projactila 

n^ 

bata 

B 

Balliatic  factor 

N/m2 

vrx 

^rx 

X  conponant  of  valocity 

Bi/sec 
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TABLE  4.2  (CONTINUED) 


SYMBOLS 

MEANING 

UNITS 

Computer 

Code 

Text 

vry 

Y  couponont  of  velocity 

m/sec 

arx 

*rx 

X  component  of  acceleration 

m/sec2 

ary 

ary 

Y  component  of  velocity 

m/sec2 

ARAD 

8 

Polar  coordinate  in  the 
orbit  plane 

Radians 

R 

R 

Radial  coordinate  in  the 
orbit  plane 

m 

<P 

Elevation  angle 

Radians 

DRIDTI 

dR/dt 

First  time  derivative  of  R 

m/sec 

DR2DT2 

d2R/dt2 

Second  time  derivative  of  R 

m/sec^ 

DAIDTI 

dA/dt 

First  time  derivative  of  A 

Rad/sec 

DA2DT2 

d2A/d2t 

Second  time  derivative  of  A 

Rad/sec 

DIP 

Ze 

Z  coordinate  of  surface 
of  earth 

km 

dqdt 

• 

q 

Rate  of  heating  per  unit 
area 

Watt/ cm 

qh 

q 

Total  heating  per  unit 
area 

J/cm^ 

erosem 

dl 

Depth  of  erosion 

cm 

TCH 

ThC 

Temperature  of  the  air 

Degrees 

Celsius 

TKH 

ThK 

Temperature  of  the  air 

Degrees 

Kelvin 

TNOSE 

To 

Temperature  at  the 
projectile  nose 

Degrees 

Kelvin 

SOS 

a 

Speed  of  sound 

m/sec 
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TABLE  6.2  (CONTINUED) 


SYMBOLS 

MEAl^ING 

UNITS 

Computer 

Code 

Teid: 

MACH 

m 

Mach  number 

Non  Dimen 

PNOSE 

Po 

Pressure  at  the  projectile 

Atm 

PH 

Ph 

Pressure  of  the  air 

Atm 

DENH 

Ph 

Density  of  the  air 

Kgm/m^ 

MU 

M 

Kinematic  viscosity  of 
the  air 

ft^/sec 

LENGTH 

L 

Typical  length 

m 

RNRB 

Rn/Pfa 

Ratio  of  radius 

Rad 

D 

D 

Aerodynamic  drag 

Newton 

S 

S 

Reference  area 

m^ 

dm 

dm 

Mass  erosion 

Kg 

ht 

ht 

Total  enthalpy 

j/Kg 

hw 

hw 

Carbon’-carbon  enthalpy 

j/Kg 

3.  Basic  Assumptions  and  Relations 
a .  General 

The  guidance  model  was  derived  from  the  reentry 
module  of  the  flight  simulation.  Most  assumptions  and 
relations  of  that  model  are  applicable  here.  Additional 
assumptions  and  relations  are  described  below. 
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b.  Additional  Assumptions  and  Relations 

The  construction  of  this  model  is  based  upon  the 
following  additional  assumptions  and  relations: 

1)  The  projectile  is  equipped  with  some  aerodynamic  means 
with  which  1  '*  ^t  can  be  produced . 

2)  The  relation  between  the  lift  coefficient  and  the  drag 
coefficient  may  be  described  in  the  form: 


Co  -  Cdo  +  KCl2 


where: 

Cq  »  The  coefficient  of  drag  of  the  projectile 

Cqq  *  The  coefficient  of  parasite  drag  (drag  at  no¬ 
lift  angle  of  attack) 

Cl  "  The  coefficient  of  lift 

K  »  The  coefficient  of  induced  drag. 

3)  The  projectile  flies  on  a  proportional  guidance  path. 

For  any  rate  of  change  in  the  orientation  of  the  line 
of  sight  between  the  projectile  and  the  target,  the 
projectile  reacts  in  a  proportional  rate  of  change  of 
its  flight  path  and  in  the  same  direction.  The 

constant  of  proportion  (Kpp)  serves  as  a  parameter  and 
is  usually  called  the  navigation  ratio. 

4)  The  lift  curve  is  symmetrical  for  a  positive  or  a 

negative  angle  of  attack. 

5)  The  projectile  responds  instantly  to  any  input  from 

the  guidance  system.  Therefore  the  projectile  always 
follows  the  desired  ideal  path  of  flight. 

6)  The  amount  of  lift  that  may  be  produced  by  the 

projectile  is  limited  by  the  stall  angle.  Therefore 
whenever  the  lift  demand  to  keep  the  desired  path  is 
beyond  the  maximum  lift  capability,  the  model  assumes 
that  the  projectile  is  flying  at  the  maximum  angle  of 
attack  but  not  beyond. 


110 


C.  THE  COMPUTER  CODE 

1.  gfflgral 

The  guidance  model  code  is  based  mainly  on  the  code 
of  the  reentry  module  of  the  flight  ^^imulation.  In  addition 
the  guidance  equations  were  modeled  and  inserted  into  the 
code. 

2.  fiov.chart 

The  flowchart  of  this  program  is  given  in  Figure 

6.1. 


D.  USING  THE  PROGRAM 

1.  InBUi; 

The  input  to  the  guidance  program  is  similar  to  the 
input  for  the  reentry  module  of  the  flight  simulation.  In 
addition  the  miss  distance  of  the  projectile  at  the  moment 
of  acquisition  of  the  target  has  to  be  entered.  Acquisition 
of  the  target  by  the  guidance  system  is  assumed  to  be 
reached  at  an  altitude  of  95-100  km  above  sea  level. 

2 .  Output 

The  output  of  the  guidance  program  is  similar  to  the 
output  of  the  launch  and  reentry  models  of  the  flight 
simulations.  The  only  exception  is  that  in  order  to  study 
the  projectile  maneuvers  the  value  of  q  is  shown  both 
graphically  (it  replaces  the  value  of  the  rate  of  heating) 
and  numerically  on  the  screen. 

A  typical  layout  of  results  from  the  guidance 
program  ic  shown  in  Figure  6.2. 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


Figure  6.2  Guidance  Layout 
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VII.  THE  COMPUTER  AND  THE  PROGRAMMING  lANGUAGES 


A.  THE  COMPUTERS 

The  mein  concept  of  this  research  was  to  create  a  handy 
and  easy  to  use  tool  for  prelixainary  design.  The  state-of- 
the-art  of  microcomputers  is  so  advanced  that  the  powerful 
CPU,  the  large  memory  and  the  excellent  graphic  capabilities 
that  most  of  them  offer  today,  seem  to  be  satisfactory  for 
this  purpose. 

A  decision  was  made  to  develop  these  models  on  the  Atari 
ST-1040  computer.  The  cost  per  performance  of  this  machine 
was  by  far  superior  to  any  alternative.  Technical 
description  of  this  machine  is  given  in  the  Appendix. 

The  new  generation  of  the  IBM  microcomputers  (lBM-2 
series)  is  similar  in  concept  to  the  Atari-ST.  These  new 
machines  use  the  3.5  inch  magnetic  discs.  The  basic 
language  of  both  machines  is  very  similar  as  well.  As  a 
result  the  transfer  of  these  models  to  the  new  IBM  machines 
is  a  relatively  short  and  easy  process. 

B.  THE  PROGRAMMING  LANGUAGE 

The  models  were  written  initially  in  ST-Basic.  This  is 
a  very  powerful  version  of  the  Basic  language.  The  big 
advantages  of  using  the  ST-Basic  language  are: 

a)  The  ST-Basic  language  is  an  interpreter-based 
language.  Programs  may  be  written  or  improved  and 
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th«n  b«  run  Inoasdiatttly  without  th«  nood  to  nako  any 
conpilation  or  linkagas. 

b)  Tha  ST-Baaic  aasuras  aaay  and  afficlant  uaa  of  tha  ST- 
aachina  and  its  graphic  capability. 

In  ordar  to  accalarato  running  tha  programs,  aftar  tha 
davalopmant  of  tha  coda  was  complatad,  a  spacial  varsion  of 
them  was  praparad  for  compilation.  Compilation  was 
complated  by  tha  LDW  compilar.  This  first  version  of  tha 
software  was  not  yet  free  from  bugs.  However,  with 
professional  assistance,  LDW  programmers  compiled  versions 
of  most  modules. 

In  order  to  ease  the  use  of  the  modules,  they  were 
packed  in  one  package.  A  managing  program,  "INTRODUl"  was 
added.  This  program  loads  and  runs  the  different  modules 
upon  the  user's  request.  A  typi..'al  layout  of  the  screen  of 
this  program  is  shown  in  Figure  7.1. 

In  order  to  save  time  and  to  prevent  the  need  to 
transfer  relevant  results  from  one  module  to  another,  the 
chaining  capabilities  of  the  software  were  used.  A 
schemati''  description  of  the  inter-modules  chaining  is  shown 
in  Figure  7.2, 

C.  RUNNING  THE  MODEL  ON  THE  MAINFRAME 

In  order  to  check  and  to  compare  the  results  of  the  new 
suggested  method  for  a  solution  of  the  "two  body  problem"  to 
the  existing  methods,  a  Fortran  version  of  the  orbit  program 
was  created.  This  program  was  run  on  the  Naval  Postgraduate 
School's  mainframe  IBM  370/3033AP  computer  in  order  to 
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Figure  7 . 1  Introduction  Layout 
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Figure  V.2  Interconnections  Between  Modules 


achittvtt  v«ry  long  runn.  Tho  advantagos  of  thm  naw  nathod 
ara  aoat  apparant  for  long  runs. 

All  othar  aodala  wara  writtan  for  aioroconputara.  Thara 
ia  no  advantaga  to  run  thaa  on  tna  aainfraaa.  In  ordar  to 
achiava  tha  graataat  banafit  froa  thaaa  aodals,  thay  hava  to 
ba  run  on  an  Atari-ST  or  on  aachinaa  with  siailar  graphics 
capabllitlas. 
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VIZI.  CA8g-STtf01BS 


GENERAL 

In  order  to  deaonatrete  the  potential  use  of  these 
sodela  for  ptt^llsinary  design,  a  fev;  case  studies  vers 
initiated.  The  following  are  the  results  of  a  few  case 
studies  based  on  these  nodels: 

1)  Parasetric  study  of  the  performance  of  a  railgun 

system 

2)  Sensitivity  of  the  range  to  errors  in  launch 
conditions 

3}  Preliminary  design  of  an  airborne  railgun  system 

4}  Preliminary  design  of  a  seaborne  railgun  system. 

B.  CASE  STUDY  1~PAJRANETRXC  STUDY  OF  THE  PERFORMANCE  OF  A 

RAILGUN  SYSTEM 

1.  Purpose 

The  aim  of  this  study  is  to  learn  about  the  effect 
of  relevant  parameters  on  the  range  of  a  projectile  fired  by 
a  railgun. 

2. 

A  nominal  case  was  selected.  In  this  example  the 
conditions  for  the  nominal  case  were  selected  as  the 
following: 

a)  Firing  of  a  50  kg  projectile 

b)  Launch  at  sea  level 

c)  The  projectile  with  a  radius  of  0.1  m 
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d)  Th«  proj«ctil«  with  a  low  drag  ahapa  (C^  *0.1) 

a)  Tha  projactlla  waa  firad  at  an  angla  of  50  dagraas 
abova  tha  local  horiion. 

Tha  flight  alaulation  aodal  waa  run  rapatitivaly. 
In  aach  run  a  diffarant  valua  of  a  aalaotad  paraaatar  waa 
ohoaan.  All  othara  wara  kapt  in  thair  noainal  valuaa.  Tha 
affaet  of  clianging  aach  paraaatar  ovar  iti?  full  ranga  ia 
ahown  by  grapha. 

3. 

Tha  raaulta  of  thia  atudy  ara  ahown  in  tha  following 


figuraa: 

a)  Pigura  8.1 — tha  affaet  of  tha  projectile  drag 
coefficient 

b)  Figure  8.3— tha  affect  of  launch  apeed 

c)  Figure  8.3— tha  effect  of  launch  altitude 

d)  Figure  8.4— tha  affect  of  launch  angla. 

4.  Concluaiona 

The  following  concluaiona  were  drawn  frost  the  study: 

a)  The  coefficient  of  drag  is  a  key  paraneter  affecting 
the  range  perfomianca.  For  axaaple,  reducing  the  drag 
coefficient  from  0.2  to  0.1  will  triple  the  range. 

b)  The  function  of  the  range  vs.  the  launch  speed  looks 
parabolic  at  the  lower  speed  range  (under  7  ksi/sec  in 
this  specific  case  study)  and  becoaes  linear  at  higher 
speeds.  This  aeans  that  the  systea  becomes  less 
efficient  as  the  launch  speed  increases. 

c)  The  range  increases  linearly  with  altitude  in  the 
najori'^y  of  tha  flight  altitudes.  The  rate  is 
approxiaately  120  NN  of  range  per  1  km  of  altitude. 
This  means  that  firing  from  an  airplane  flying  at 
typical  jet  cruise  altitude  of  about  36,000  ft  may 
double  the  range  in  comparison  with  the  range  achieved 
by  firing  at  sea  level. 
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RePROOUCED  AT  GOVERNMENT  EXPENSE 


EFFECT  OF  DRAG  COFFICIENT 


Figure  8 . 1 


The  Effect  of  the  Projectile  Drag  Coefficient 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


EFFECT  OF  LAUNCH  SPEED 


e 


Figure  8.2  The  Effect  of  Launch  Speed 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


EFFECT  OF  ALTITUDE 
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Figure  8.3  The  Effect  of  Launch  Altitude 
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REPRODU<;ED  AT  GOVERNMENT  EXPENSE 


EFFECT  OF  LAUNCH  ANGLE 
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Figure  8.4  The  Effect  of  Launch  Angle 
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d)  From  th«  graph  of  the  range  vs.  launch  angle,  it  can 
be  learned  that  firing  all  high  angles  (above  45^) 
gives  longer  ranges.  This  is  caused  by  the  effect  of 
the  aerodynamic  drag  in  the  atmosphere.  Firing  at  a 
higher  angle  results  in  lengthening  the  phase  of  space 
flight  and  shortening  the  phase  of  atmosphere  flight. 
This  is  a  more  efficient  use  of  energy — less  kinetic 
energy  is  lost  to  the  atmosphere.  In  this  case  study, 
maximum  range  was  achieved  at  an  angle  of  launch  of 
about  52  degrees. 


C.  CASE  STUDY  2— SENSITIVITY  OF  THE  RANGE  TO  ERRORS  IN 

LAUNCH  CONDITIONS 

1.  Aim 

The  aim  of  this  study  was  to  check  the  sensitivity 
of  the  range  to  errors  in  launch  parameters.  In  addition, 
the  demands  from  the  guidance  system  to  correct  such  errors 
were  exeunlned. 

2.  Method 

In  order  to  study  the  sensitivity  of  the  range  to 
errors  in  launch  parameters,  the  flight  simulation  code  was 
run  for  the  following  cases: 

a)  Launch  at  the  nominal  condition 

b)  Launch  with  a  small  deviation  from  the  nominal  launch 
angle 

c)  Launch  with  a  small  deviation  from  the  nominal  launch 
speed 

d)  Launch  of  a  projectile  that  differs  slightly  in  its 
coefficient  of  drag  from  the  nominal  value. 

In  order  to  study  the  demands  from  the  guidance 
system  it  was  asstimed  that  at  reentry  the  guidance  system 
discovers  a  range  error.  The  measure  of  the  range  error  was 
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chosen  as  the  sum  ot  all  the  errors  caused  by  the  deviations 
from  the  nominal  launch  conditions. 

The  guidance  system  controlling  the  flight  path  of 
the  projectile  corrects  the  range  errors.  The  output  from 
the  guidance  program  includes  a  graph  of  the  lift  forces 
needed  to  achieve  the  desired  corrections.  The  demand  for 
lift  may  be  used  to  design  the  lift  surfaces  and  the 
guidance  system  of  the  projectile. 

3.  RfiAttita 

The  results  of  this  study  are  summarized  in  Table 
8.1.  In  this  table  the  sensitivity  coefficients  of  the 
range  to  errors  in  launch  conditions  and  to  deviation  of  the 
projectile  from  its  design  specification  are  presented. 
These  coefficients  were  calculated  for  the  following  nominal 
case: 

a)  Projectile  mass  (Mp)~50  kg 

b)  Coefficient  of  drag  (C^j)— 0.096 


TABLE  8.1 

SENSITIVITY  COEFFICIENTS 


The  Cause 

Deviation  of  1  m/s  in  launch  speed 

Deviation  of  1  degree  in  the  angle 
of  launch 

Deviation  of  0.01  in  the  drag 
coefficient 


Range  Error 
9.15  NM 
5.80  NM 

8.09  NM 
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c)  Launch  angle — 50® 

d)  Muzzle  speed— 7  kin/s. 

To  correct  such  errors  the  projectile  must  have  a 
turn  capability  of  7.5  g's.  This  Is  the  result  of  a  run 
made  with  the  guidance  model .  The  run  was  made  for  a 
proportional  navigation  guidance  system  with  coefficient  of 
proportion  of  3.0.  Greater  refinement  In  the  production  of 
the  projectile  and  less  deviation  In  the  launch  conditions 
may  ease  the  demand  from  the  guidance  system. 

4.  Conclusions 

In  order  to  achieve  greater  effectiveness  of  the 
rallgun  weapon  system,  the  warhead  has  to  be  brought  as 
close  as  possible  to  the  target.  .he  measure  of  accuracy 
needed  is  dependent  on  the  type  of  warhead  and  Its  yield. 

The  range  Is  very  sensitive  to  errors  In  launch 
conditions  and  to  Inaccuracies  In  the  projectile 
manufacturing.  Increasing  the  accuracy  of  the  launching 
system  is  not  cost  effective. 

Whenever  greater  accuracy  Is  desired  a  guidance 
system  Is  needed.  A  proportional  guidance  system  may 
correct  these  errors  and  assure  high  effectiveness  of  the 
overall  performance  of  the  rallgun  system. 


D.  CASE  STUDY  3— SEABORNE  SYSTEM 


1.  Purpoaa 

Thtt  aim  of  this  study  was  to  check  the  feasibility 
and  to  estimate  the  performance  of  a  railgun  system  mounted 
aboard  a  naval  vessel. 

2. 

By  using  the  flight  simulation  program  the  basic 
parameters  of  the  railgun  system  were  defined  in  order  to 
assure  the  desired  ranges  and  effectiveness. 

Once  the  basic  parameters  (mass  and  coefficient  of 
drag)  were  defined,  the  basic  shape  of  the  projectile  was 
identified  by  using  the  drag  estimation  module. 

The  results  of  the  drag  module  and  the  flight 
simulation  were  fed  into  the  railgun  model.  The  output  of 
the  railgun  program  included  the  cross  section  of  the  barrel 
shape  and  estimates  of  the  overall  system  weight  and  volume. 

3 .  Results 

The  results  of  this  study  are  presented  in  the 
following  figures: 

a)  Figure  8.5— 'final  output  of  the  flight  simulation 
model 

b)  Figure  8.6—the  warhead  shape  derived  by  the  drag 
module 

c)  Figure  8.7 — results  of  the  railgun  simulation 

d)  Figure  8.8 — the  railgun  structure 

e)  Figure  8.9~summary  of  the  power  conditioning  system 

f)  Figure  8.10 — total  weight  of  the  cooling  system  and 
its  main  subsystems 
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Figure  8.6  Case  Study  3— The  Warhead  Shape 
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Figure  8.7  Case  Study  3— The  Railgun  Simulation 
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Figure  8.8 


RAIL-SQI  STRUCTQRE 
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Figure  8.9  Case  Study  3 — The  Power  Conditioning 
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Figure  8.10  Case  Study  3 — The  Cooling  System 


q)  Figure  8.11— >auHi«ry  of  tho  woight  and  alia  of  tha 
whola  alrborna  systaa. 

4.  consluiioni 

Tha  aaln  concluaion  tram  this  praliainary 
axaaination  is  that  it  is  faasibla  to  dasign  a  railgun 
systaa  for  a  naval  vassal.  Such  a  systaa  %rould  ba  eapabla 
of  firing  a  SO  leg  projactila  to  a  ranga  on  tha  ordar  of  lOOO 
NN.  This  parforaanea  is  aasurad  if  tha  valua  of  tha  drag 
coafficiant  is  0.1  or  lass.  This  low  valua  can  ba  achiavad 
by  propar  dasign  of  tha  projactila  aarodynaaic  shapa. 

Tha  railgun  dasign  will  ba  basad  on  a  20  ca  (width) 
rail.  Tha  coppar  rails  will  hava  a  2  by  6  ca  rectangular 
cross  section.  Tha  railgun  length  has  to  ba  at  laaat  8.0  a 
long. 

Tha  inductive  power  conditioning  systaa  will  weigh 
about  1116  tons  and  its  voluaa  will  ba  1610  cubic  asters 
(Figure  8.9) . 

E.  CASE  STUDY  4— AIRBORNE  SYSTEM 

1.  purpoat 

The  ala  of  this  study  was  to  axaalna  the  feasibility 
and  the  perforaance  of  an  airborne  railgun  systea  and  to 
obtain  an  estlaate  of  its  weight  and  size  in  order  to  check 
the  capability  of  airplanes  to  carry  tha  systcun. 

2. 

The  aethod  of  this  exaaination  was  very  siailar  to 
the  exaaination  of  the  naval  railgun  systea  (Case  3) .  In 


135 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


1 


Desk  File  Run  Edit  Debus 


Figure  8,11  Case  Study  3 — Total  Weight  and  Size 
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this  case  typical  airborne  parameters  were  used  by  the 
airborne  modules. 

3.  Results 

The  results  of  this  study  are  presented  in  the 
following  figures: 

a)  Figure  8.12 — final  output  of  the  flight  simulation 
model 

b)  Figure  8.13-- the  warhead  shape  derived  by  the  drag 
module 

c)  Figure  8.14 — results  of  the  railgun  simulation 

d)  Figure  8.15 — the  railgun  structure 

e)  Figure  8.16 — summary  of  the  power  conditioning  system 

f)  Figure  8.17 — total  weight  of  the  cooling  system  and 
its  main  subsystems 

g)  Figure  8. IS— summary  of  the  weight  and  size  of  the 
whole  airborne  system. 

4.  sposl-ttaion 

The  conclusion  of  this  preliminary  study  is  that  it 
appears  feasible  to  install  a  railgun  system  on  an  existing 
airplane.  This  system  would  be  capable  of  firing  a  50  kg 
projectile  to  a  range  on  the  order  of  200  nautical  miles 
(assuming  that  a  C(j[  of  about  0.1  will  be  achieved  by  proper 
design  of  the  aerodynamic  shape  of  the  projectile) . 

The  railgun  system  will  be  based  on  2*16  cm 
rectangular  copper  rails  with  a  span  of  about  20  cm.  The 
inductive  power  system  weights  about  70  tons  and  its  volume 
is  about  37  cubic  meters. 
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Figure  8.14  Case  Study  4 — Railgun  Simulation 
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IHiCTiyE  POWER  COHDiTIOHiRS  SVSTER 

THE  HElSrorTHriiSoPOLAR-KiiERftTOr’ir  47 . 151  TON 

THE  yOLlftE  OF  THE  HOnOPOLAR-SEHERATOR  IS  23.576  H''3 

THE  WEIGHT  OF  THE  INDUCTIUE  SVSTEII  IS  57.526  TON 

THE  miHE  OF  THE  INDOCTIUE  SVSTEII  IS  113.482  ;t^3 

THE  WEIGHT  OF  THE  SWITCHING  SVSTEII  IS  4.584  TON 

THE  vOLlUE  OF  THE  SWITCHIHS  SVSTEII  IS  3.751 
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Figure  8.16  Case  Study  4 — Power  Weight  and  Size 
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Figure  8.17  Case  Study  4 — Cooling  Weight  and  Size 
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Figure  8.18  Case  Study  4 — Total  Weight  and  Size 
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Thtt  overall  size  and  weight  of  the  system  is  about 
90  tons.  It  appears  that  a  railgun  system  of  this  weight 
and  size  is  suitable  to  be  mounted  on  an  existing  airplane. 
This  specific  system  will  require  an  airplane  the  size  and 
performance  of  a  Boeing  747.  Longer  ranges  may  be  achieved 
if  the  system  is  installed  on  larger  planes  in  the  future. 

Firing  the  projectile  at  a  rate  of  10  projectiles 
per  minute  will  demand  a  power  supply  on  the  order  of  33.5 
Mwatts.  The  wasted  heat  will  be  cooled  by  large  radiators 
mounted  on  top  of  the  aircraft.  An  imaginary  picture  of 
such  an  airplane  is  shown  in  Figure  £.19. 


IX.  RXCQMMENDXTIQWS 


My  attitud*  in  craating  this  nodal  was  to  produca  an 
afficiant  tool  for  pralininary  dasign.  For  tha  purposa  of 
pralininary  dasign  fast  raaction  tins  is  mora  iaportant  than 
graatar  accuracy,  othar  nodals  do  axist.  sona  of  than  ara 
very  sophist icatad  and  taka  aora  tiaa  to  run.  It  is 
raconmandad  that  a  conparison  of  tha  rasults  of  this  nodal 
with  tha  rasults  of  othar  nodals  be  parfomad. 

Tha  fiald  of  railgun  and  its  auxiliary  aystans  is 
advancing  quickly.  Many  efforts  ara  nada  to  inprova  tha 
railgun  parfomanca  and  its  technology.  Tha  intensive 
research  and  tha  variolas  projects  in  progress  nay  produce 
new  infomation  and  data.  It  is  raconnandad  that  new 
infomation  and  updated  data  be  continuously  inplenanted 
into  tha  nodal. 

It  is  raconnandad  that  additional  nodules  vill  be  added 
to  this  package.  The  following  topics  have  to  be  covered: 

a .  Weapons : 

1)  High  energy  lasers 

2)  Particle  bean 

3)  Kinetic  energy  rockets 

b.  Battle  nanagenent 

1)  Sensors 

2 )  Tracking 
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3)  Coaputlng 

4)  Conunication 
e.  Non-soz  alsaions 

1)  Anti-aircraft 

2)  Anti-aurfaca. 
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